
m M 

: yyf. 


m : 


; . ■ r."' ; V. : y 


MSI i 


■ ! « 

•'-■•• . ; 


BB wis 


r <>mWk‘ " 

Bill 

Wiii* 


: ., 

■ ■ ; • :-' ■■•■>• 

’ ■ • • • : : • 

B BB 4 BB 


iiil! 

* \ 




'■A M 

■ ;•;!>: : 

: ; 


MRNI 






b 4 f 


• .': ; 'B -;i 


V> ' ' - : . 


; w1 


r 


ASHIN6T0M, I 


1. Report No. 2. Government Accession No. 

NASA CR-2548 

3. Recipient's Catalog No. 

4. Title and Subtitle 

STEADY -STATE AND TRANSIENT ANALYSIS OF A SQUEEZE 
FILM DAMPER BEARING FOR ROTOR STABILITY 

5. Report Date 

May 1975 

6. Performing Organization Code 

7. Author(s) 

L. E. Barrett and E. J. Gunter 

8. Performing Organization Report No. 

ME -4040-114 -73U 

10. Work Unit No. 

9. Performing Organization Name and Address 

University of Virginia 
Charlottesville, Virginia 22902 

1 1 . Contract or Grant No. 

NGL 47-005-050 

13. Type of Report and Period Covered 

Contractor Report 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


Final Report. Project Manager, Robert E . Cunningham, Fluid System Components Division, 
NASA Lewis Research Center, Cleveland, Ohio 


16. Abstract 

This report presents a study of the steady-state and transient response of the squeeze film 
damper bearing. Both the steady -state and transient equations for the hydrodynamic bearing 
forces are derived. The steady-state equations are used to determine the bearing equivalent 
stiffness and damping coefficients. These coefficients are used to find the bearing configura- 
tion which will provide the optimum support characteristics based on a stability analysis of the 
rotor -bearing system. The effects of end seals and cavitated fluid film are included. The 
transient analysis of rotor -bearing systems is performed by coupling the bearing and journal 
equations and integrating forward in time. The effects of unbalance, cavitation, and retainer 
springs are included in the analysis. Methods of determining the stability of a rotor -bearing 
system under the influence of aerodynamic forces and internal shaft friction are discussed. 
Particular emphasis is placed on solving the system characteristic frequency equation, and 
stability maps produced by using this method are presented. The study shows that for optimum 
stability and low force transmissability the squeeze bearing should operate at an eccentricity 
ratio e < 0. 4. 


17. Key Words (Suggested by Author(s)) 

Damper bearing; Oil squeeze film; Steady- 
state analysis; Transient analysis; Rotor 
stability; Rotor dynamics 

18. Distribution Statement 

Unclassified - unlimited 
STAR category 07 (rev. ) 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of Pages 

107 

22. Price* 

$5.25 


For sale by the National Technical Information Service, Springfield, Virginia 22151 
















TABLE OF' CONTENTS 


CHAPTER 1 INTRODUCTION 1 

CHAPTER 2 THEORETICAL ANALYSIS 

2.1 Reynolds Equation ..... 4 

2.2 Bearing Forces in Fixed Coordinates . 6 

2.3 Bearing Cavitation 10 

2.4 Bearing Forces in Rotating Coordinates 12 

CHAPTER 3 ROTOR-BEARING STABILITY AND STEADY- STATE 

ANALYSIS 

3.1 Rotor-Bearing Stability ....... 23 

3.2 Steady-State Analysis 33 

CHAPTER 4 TRANSIENT ANALYSIS 

4.1 Introduction 42 

4.2 Analysis ........ 44 

CHAPTER 5 CONCLUSIONS AND RECOMMENDATIONS 

5.1 Predicting Rotor Instability ..... 71 

5.2 Determining the Stiffness and Damping 
Coefficients of the Squeeze Film 

Damper Bearing ............ 72 

5.3 Transient Analysis .......... 73 

5.4 Advantages of Bearing Simulation . . . 74 

5.5 Limitations of Analytical Investigations 74 

5.6 Recommendations for Future Research . . 75 

BIBLIOGRAPHY 77 

APPENDIX A DESCRIPTION OF PROGRAM SQFDAMP 79 

iii 



LIST OF FIGURES 


FIGURE PAGE 

2-1 Squeeze Film Damper Bearing Configuration in 5 

Fixed and Rotating Coordinate Systems 

2-2 Comparison of Finite Length and Short Bearing 7 

Solutions 

2-3 Uncavitated Pressure Profile Showing Region 11 

of Negative Hydrodynamic Pressure 

2-4 Cavitated Pressure Profile 'With Negative 13 

Hydrodynamic Pressure Equated to Zero 


2-5 k. Axial Pressure Distribution of Bearing With 19 

Circumferential Oil Groove 

b. Axial Pressure Distribution of Bearing With 
End Seals and Circumferential Oil Groove 


3-1 Three-Mass Flexible Rotor Mounted On Flexible, 26 

Damped Supports 

3-2 Characteristic Matrix For Three-Mass Model In- 27 

eluding Aerodynamic Cross-Coupling, Internal 
Shaft Friction and Absolute Shaft Damping 

3-3 Stability of a Flexible Rotor With Aerodynamic 30 

Cross Coupling (Q = 20,000 lb /in, N - 10,000 RPM) 

3-4 Stability of a Flexible Rotor With Aerodynamic 31 

Cross Coupling (Q = 100,000 lb/ ini, N =10,000 RPM) 

3-5 Damping Coefficient For Squeeze Film Bearing 35 

With Cavitated Film - End Seals And Oil Supply 
Groove Included 


3-6 Stiffness Coefficient For Squeeze Film Bearing 36 

With Cavitated Film - End Seals And Oil Supply 
Groove Included 

3-7 Maximum Pressure For Squeeze Film Bearing With 37 

Cavitated Film - End Seals And Oil Supply 
Groove Included 

3-8 Damping Coefficient For Squeeze Film Bearing 39 

With Uncavitated Film - End Seals And Oil 
Supply Groove Included 


iv 



FIGURE PAGE 

3- 9 Maximum Pressure For Squeeze Film Bearing 40 

With Uncavitated Film - End Seals And Oil 
Supply Groove Included 

4- 1 Unbalanced Rotor In Cavitated Squeeze Film 46 

Bearing - L = 0,45 IN. - Unbalance Eccentricity 
= 0.002 IN. - No Retainer Spring 

4-2 Stiffness Coefficient For Squeeze Film Bearing 49 

Of Figure (4-1) 

4-3 Damping Coefficient For Squeeze Film Bearing 50 

Of Figure (4-1) 

4-4 Unbalanced Rotor In Cavitated Squeeze Film Bear- 51 

ing L = 0.90 IN. - Unbalance Eccentricity = 

0.002 IN. - No Retainer Spring 

4-5 Stiffness Coefficient For Squeeze Film Bearing 52 

Of Figure (4-4) 

4-6 Damping Coefficient For Squeeze Film Bearing Of 53 

Figure (4-4) 

4-7 Unbalanced Rotor In Cavitated Squeeze Film Bear- 55 


ing - L = 0.45 IN. - Unbalance Eccentricity = 

0.002 IN. - Retainer Spring Stiffness, KR - 
123,000 LB/IN. 

4-8 Unbalanced Rotor In Cavitated Squeeze Film Bear- 56 

ing - L = 0.90 IN. - Unbalance Eccentricity = 

0.001 IN. - No Retainer Spring 

4-9 Unbalanced Rotor In Cavitated Squeeze Film Bear- 57 

ing - L = 0.90 IN - Unbalance Eccentricity = 

0.001 IN. - Retainer Spring Stiffness, KR = 

123,000 LB/IN 

4-10 Unbalanced Rotor In Cavitated Squeeze Film Bear- 59 

ing - L = 0.90 IN. - Unbalance Eccentricity - 
0.002 IN. - Retainer Spring Stiffness, KR = 

123,000 LB/IN 

4-11 Unbalanced Rotor in Uncavitated Squeeze Film Bear- 60 

ing - L = 0.90 IN. - Unbalance Eccentricity = 

0.002 IN. - Retainer Spring Stiffness, KR = 

123,000 LB/IN 


v 



FIGURE 


PAGE 


4-12 Vertical Unbalanced Rotor In Squeeze Film Bear- 61 

ing - Effect Of Unbalance Magnitude - Unbalance 
Eccentricity - 1.75 Mils 

4-13 Vertical Unbalanced Rotor In Squeeze Film Bear- 62 

ing - Effect Of Unbalance Magnitude - Unbalance 
Eccentricity - 2.10 Mils 

4-14 Vertical Unbalanced Rotor In Squeeze Film Bear- 63 

ing - Effect of Unbalance Magnitude - Unbalance 
Eccentricity = 2,45 Mils 

4-15 Vertical Unbalanced Rotor In Squeeze Film Bear- 64 

ing - Effect of Unbalance Magnitude - Unbalance 
Eccentricity = 3,50 Mils 

4-16 Horizontal Unbalanced Rotor In Squeeze Film Bear- 66 

ing - Effect Of Retainer Springs - Retainer 
Spring Stiffness, KR - 0 LB/IN 

4-17 Horizontal Unbalanced Rotor In Squeeze Film Bear- 67 

ing - Effect Of Retainer Springs - Retainer 
Spring Stiffness, KR = 50,000 LB/IN 

4-18 Horizontal Unbalanced Rotor In Squeeze Film Bear- 68 

ing - Effect Of Retainer Springs - Retainer 
Spring Stiffness, KR = 100,000 LB /IN 

4-19 Horizontal Unbalanced Rotor In Squeeze Film Bear- 69 

ing - Effect of Retainer Springs - Retainer 
Spring Stiffness, KR = 200,000 LB /IN 


vi 



NOMENCLATURE 


SYMBOL 

C 

Co 

c , c 

xx xy 
C , C 

yy yx 

e 

EMU 


Fx 

Fy 

F r 

F e 

FMAX 

FU 


FURATIO 


h 



k 

o 


KR 




N 



DESCRIPTION 
Bearing Clearance 
Equivalent bearing damping 


UNITS 

in. 


lb-sec /in 


Bearing damping Ib-sec/in 


Journal Eccentricity in. 

Ratio of unbalance eccentricity 
to bearing clearance 

Force component in x-direction lbf 

Force component in y-direction lbf 

Force component in radial direction lbf 

Force component in tangential direction lbf 

Maximum hydrodynamic force lbf 

Force due to rotating unbalance lbf 

Ratio of FMAX to FU - — - 

Fluid film thickness in. 

Unit vectors in fixed coordinate system 

Equivalent bearing stiffness lb /in 

Retainer Spring Stiffness lb /in 


Bearing Stiffness 


lb / in 


Bearing length in. 

Rotor speed RPM 

Unit vecotors in rotating coordinate 

system 


vii 



SYMBOL 

DESCRIPTION 

UNITS 

P 

Pressure 

lb/in2 

PMAX 

Maximum hydrodynamic pressure 

lb/in^ 

Q 

Aerodynamic cross coupling 

lb /in 

R 

Bearing radius 

in 

t 

Time 

sec 

W 

Weight 

Ibf 

x,y,z 

Displacements 

in. 

y 

Viscosity 

lb-sec 

e,e' 

Angular measure 

— 

1 

Journal precession rate 

sec"-*- 

or 

Angular velocity 

sec~l 


xri i i 



CHAPTER 1 


INTRODUCTION 

Modern turbomachines are highly complex systems. Current 
design trends are producing machines that consist of several 
process stages joined together. The rotors in these machines 
are highly flexible shafts, often mounted in more than two bear- 
ings, that rotate at very high speeds. It is not uncommon to 
see machines that operate above the second critical speed. As 
a result the system dynamics are very complicated. 

One of the major problems encountered in these machines is 

instability produced by aerodynamic forces on impeller wheels, 

friction in the stressed rotor and hydrodynamic forces in the 
bearings. The instability is characterized by large amplitude 

whirl orbits and often results in bearing or total machine failure. 

It is often aggravated by unbalance and other external forces 

transmitted to the machine. Production losses from failed machines 

are very high and it may take many months to repair or replace the 

failed unit. In addition operator safety is jeopardized when 

machines fail and occasional loss of life occurs. 

From the earliest investigations of rotor instability, it has 
been known that the use of flexible, damped supports has an effect 
on instability and can eliminate it or alter the speed at which 
it occurs. Recent research has produced a large body of knowledge 
on the use of these supports and their effect on instability. 

The squeeze film damper bearing is one type of flexible 
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support that is currently being investigated. This study ex- 
amines the squeeze bearing and through computer simulation shows 
its effects on several rotor-bearing systems. The equations for 
the hydrodynamic bearing forces are developed in both fixed and 
rotating coordinate systems. The use of two coordinate systems 
allows for both steady-state and transient analysis of bearing 
performance. This results in more efficient bearing analysis 
and a savings in time and money when experimental testing of the 
bearings is conducted. 

The steady-state behavior of the bearing results in the for- 
mulation of bearing stiffness and damping coefficients which can 
be used to set the bearing configuration. This is accomplished 
by comparing the coefficients with required values obtained from 
a stability analysis of the rotor-bearing system. Several methods 
of determining the system stability are discussed. The effects of 
end seals and cavitation of the fluid film are also included in 
the steady-state coefficients. 

The transient analysis is very useful in determining the 
bearing response to particular forms of external and internal 
forces as noted previously. Also the effect of bearing retainer 
springs and fluid film cavitation can be found. The transient re- 
sponse is found by tracking the journal motion forward in time 
by integrating the equations of motion under the influence of the 
system forces. 

The limitations of and assumptions used in deriving the 
steady-state and transient equations are discussed in order to 
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obtain meaningful interpretation of the results and to es- 
tablish useful design criteria. 

Dr. R. Gordon Kirk developed the computer programs used 
to perform the transient analysis in this report. 
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CHAPTER 2 

THEORETICAL ANALYSIS 
2.1 REYNOLDS EQUATION 

The configuration of the squeeze film damper bearing 
is shown in Figure (2-1) where the clearance has been exag- 
gerated. Both fixed and rotating coordinate systems are shown, 
and the bearing equations are derived for both systems. The 
definitions of the various parameters are listed in the nomen- 
clature section of this report. 

The basic bearing equation is the Reynolds equation which 
is derived from the Navier-Stokes equations for incompressible 
flow. With the proper bearing parameters the equation for the 
fluid film forces are derived. [1] 

The Reynolds equation for the short, plain journal bearing 
is given in both fixed and rotating coordinates by: 

Fixed coordinates: 


L_ fill 1Z 

3Z [6y 3Z 


(a> b + uj) 


3h , 2 3h 

30 at 


Rotating coordinates: 


3 

3Z 



(to, + co. - 2^) 

■ -J 


9h , 0 3h 

TeF ^ at 


(2-1) 


(2-2) 


As shown in Figure (2-1), the angle 9 in the fixed coor- 
dinate expression is measured from the positive x-axis in the 
direction of rotation whereas the angle 9 ' in the rotating co- 
ordinate expression is measured from the line of centers in the 
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direction of rotation. The assumptions used in the derivation 
of equations (1) and (2) include: 

1. The fluid inertia terms in the Navier-Stokes equations 
have been neglected due to their small magnitude. 

2. Body forces in the fluid film have been neglected. 

3. The fluid viscosity is constant. 

4. The flow in the radial direction has been neglected, 
that is, the short bearing approximation has been used. 

Figure (2.2) shows a comparison of the short bearing sol- 
ution and the general solution of the Reynolds equation solved 
by a finite difference technique for the plain journal bearing 
under steady state conditions. It can be seen that the short 
bearing solution is highly accurate for a wide range of eccent- 
ricities for L/D < 1/4 and is acceptable for L/D values up to 1 
if the eccentricy ratio is low. The normal design range of the 
squeeze film bearings will be L/D < 1/2 and eccentricity ratios 
< 0.4. 

By assuming the bearing is perfectly aligned (h not a 
function of Z) equations (1) and (2) are integrated directly to 
yield expressions for the fluid film forces. 

2.2 BEARING FORCES IN FIXED COORDINATES 

For the plain bearing with full end leakage the appropriate 
boundary conditions are: 

(2-3) 

P(6,0) = P(0,L) = 0 
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In the fixed coordinate system the film thickness, h, is given 
by: 


h = c - x cos 0 - y sin 0 


Substituting into equation (2-1) and integrating yields: 


P (0 , Z) = — y Z 2 - LZ 


z» - Lzj + “j 


) ih + 2 £h 

30 3t 


Differentiation of equation (2-4) yields: 


36 - x sin 0 - y cos 0 


x cos 0 - y sin 0 


The incremental force acting on the journal is: 


tF = -P ( 0 , Z ) Rd0dZ \n 


The relationship between the unit vectors in the fixed and 
rotating reference frames is: 


i = cos 0 ]n r - sin 0 p Q 
J * sin 0 (n r + cos 0 |n 0 


n r - cos 0 i + sin 0 j 


n a = -sin 0 i + cos 0 


( 2 - 12 ) 
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Substituting equation (2-11) into equation (2-8) yields: 


AF = - P (0 , Z) Rd0dZ (cos 0 l - sin 0 


j) 


(2-13) 


Thti elemental x and y force components may be found by 
taking the dot product of equation (2-13) with unit vectors in 
the x and y directions. 


AF' = (AF • 1)1 = - (P (0 , Z) RdQdZ cos 0) i 

X 


(2-14) 


AFy = (AF • 3)3 


-(P(0,2)RdedZ sin 0) 3 


(2-15) 


The total force components in the x and y directions are found by 
integrating equations (2-14) and (2-15) over the entire journal 
surface. 


F 


x 



L 

P (0 , Z ) R cos 0 dZ d0 


(2-16) 


27T l 

■ - // p(9 ' 
0 0 


Z)R sin 0 dZ d0 


(2-17) 


9lx 

Substituting the expressions for — and ~ into the pressure 

equation and integrating around the bearing circumference gives: 
2tt 


_ -yRL 3 / (xsin0-ycos0) -2 (xcos0+ysin0) 


(c - x cos 0 - y sin 0) 3 


cos 0 | d 0 (2-18) 
sin 0 I 


9 



The above equation is applicable to the evaluation of 
the forces developed in the plain journal bearing as well as 
the squeeze film damper bearing for arbitrary values of journal 
dispacement, velocity, and shaft and bearing housing angular 
velocities. Hence the analysis can also be used for the general 
floating bush bearing with rotation. 

For the case of the squeeze film damper where the journal and 
housing are constrained from rotating, ( 0 ^= 0 ^ = 0), the force 
expressions become 


F xl 

fl 

— up't. 3 f -2 (X. cos 0 + y sin 0 ) 

COS 0j 

( F H 

2 J (c - x cos 0 - y sin 0) 3 
0 

^sin 0J 


d0 


(2-19) 


These non-linear fluid film forces are easily combined with 
the rotor-bearing system dynamical equations providing a complete 
non-linear dynamical analysis of the system. Because the bear- 
ing force equations are written in fixed Cartesian coordinates 
a transformation from one coordinate system to another is not 
required. This is very important for conservation of computation 
time since the bearing pressure profile must be integrated at 
each time step of the system motion. 


2.3 BEARING CAVITATION 


If the complete pressure profile is calculated without re- 
gard to cavitation or rupture of the film, then the bearing pressure 
will be similar to Figure (2.3). This figure represents the three 
dimensional pressure generated in the bearing. 

The exact mechanism causing cavitation in fluids is not 
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fully known. It is known that film rupture is influenced by gas 
and solid content of the fluid. Recent investigations have 
shown that a fluid may stand large tensile stresses [2]* 
and its ability to withstand rupture is dependent on its past 
history. In this investigation it is assumed that cavitation 
occurs when the pressure in the film drops below ambient pressure. 
The eavitated film then extends over only a section of the bear- 
ing circumference as shown in Figure (2-4). Recent experimental 
research has shown that cavitation in the squeeze bearing occurs 
in streamers of bubbles which extend around the entire bearing 
[3]. These streamers initially appear at the center of the bear- 
ing and extend outward as the rotor speed increases. It is be- 
yond the scope of this present research to analyse this type of 
cavitation effect. Therefore the conventional eavitated film 

is assumed to occur when P<P where P is the assumed cavitation 

c c 

pressure. 

When evaluating the integral of equation (2-19) , negative 
pressures are equated to zero if the film is assumed to cavitate. 
If the oil supply pressure is sufficiently high and suitable 
operating conditions exist the film does not cavitate. 

2.4 BEARING FORCES IN ROTATING COORDINATES 

The Reynolds equation in rotating coordinates was given 
by equation (2-2) . Assuming steady-state circular synchronous 
precession of the journal about the bearing center and no axial 
misalignment, equation (2-2) can be integrated in closed form. 

The resulting equations for the bearing forces give the equivalent 
stiffness and damping of the bearing. 
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Figure 2-4 Cavitated Pressure Profile With Negative 
Hydrodynamic Pressure Equated to Zero 




Applying the boundary conditions: 


p(e,o) - p(e,L) = o 


(2-20) 


and integrating equation (2-2) yields: 


p(e,z) = 3E 
h 3 



(2-21) 


where 6' is measured from the line of centers in the direction of 
rotation as shown in Figure (2-1), 

In rotating coordinates the film thickness, h, is given by: 


h 


+ e cos 0 ') 


( 2 - 22 ) 


where the eccentricty ratio, e, is defined as: 

£ = « 
c 

Differentiation of equation (2-22) gives: 
0h 

30-1= -c e sin 0' 

and 

3h 

“ ce cos 0 * 


(2-23) 


(2-24) 


(2-25) 


Substituting the expressions for h and its 
equation (2-21) and integrating gives the fluid 

V 

0 2 

F = PRL 3 j OP e sin 0 1 + e cos O') 

c 2 J (1 + e cos 9') 3 

f 

6i 


derivitives into 
film force: 


<39 ' |n r (2-26) 
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(2-27) 


The transformation into the ri 0 coordinates is 

|n = -cos 0‘ n -sin 0 ' n Q 
r r 0 

The components of the fluid film force in the radial and 
tangential directions, n r and ng are found by taking the dot 
product of the force with unit vectors in the radial and tan- 
gential directions. Thus: 

- + + + (2-28) 

F r * < p * n r )n r 


and 



n 0^ n 0 


(2-29) 


and the force components are: 

i 



/•02 




F J 

F n 1 

. = -yRL 3 | (.<j> g. sin 0 ' 

+ e cos 8 ' ) 

Jcos 

0 ' 1 d0 ' (2-30) 

l ej 

C 2 1 / -i , 

1 (1 + £ COS 

CD 

i 

[sin 

ej 


^ el 





The limits of integration, 9f and 02 , define the area over 


which a positive pressure profile exists and are dependent on 
the type of journal motion and whether or not cavitation occurs. 

It is assumed that the journal is precessing in steady- 
state circular motion about the origin. Therefore e = 0 and 
the pressure expression, equation ( 2 - 21 ) becomes: 


P ( 6 * ,Z) 


6yuoe sin 8' 

c 2 (1 + e cos 0') 3 




(2-31) 
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The maximum pressure in the axial direction occurs at Z = L/2 
and equation (2-31) is rewritten as: 


P(0',L/2) - si » 6' 

2c 2 (1 + e cos 0 ' ) 3 


(2-32) 


By differentiating equation (2-32) with respect to 0 ' and 
equating to zero, the tangential location of the maximum pressure 
may be found. The angle, 8' max., where the pressure is maximum 
is given by: 


(1 + e cos 6 ' cos 6' + 3 e sin* 6* = 0 

max max D max 


(2-33) 


The angle 6^ ax varies with e and shifts from 


0 * 

max 


3ir 


when e = 0 


to 


0 1 * ir when e = 1 

max 

and the maximum pressure is given by: 


-3yL 2 o)e sin 0 1 


max 


max 


2c 2 (1 + e cos 0 ' ) 3 

max 


(2-34) 


The pressure expressed by equation (2-32) is positive over 
the region 8' = ir to 8’ = 2ir and for the cavitated fluid film 
the limits of integration in equation (2-30) are: 


0 1 

1 


IT f 0 ' = 2tt 

2 


(2-35) 
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The radial and tangential components of the fluid film force 
are given by: 


2n 


9 


_ -yRL 3 eq) 

/~t 2 


/ sin e* 
(1 + e < 


cos e 1 


de 


(2-36) 


cos 6 ' ) 3 sin 6 

**r \ 

The integrals in equation (2-36) were integrated using Booker's 
method [4]. The resulting force components are: 


_ _ -2yRL 3 etoe (2-37) 

r c 3 ( 1 - e 2 ) 2 

and 

F e = -yRL 3 Treo) (2-38) 

2c 3 (1 - e 2 ) »/ 2 

The force in equation (2-37) appears as a stiffness coefficient 
times a displacement acting in line of the displacement towards 
the. bearing center. The equivalent bearing stiffness is: 

v _ 2yRL 3 eo) (2-39) 

o — ; — 

c 3 (1 - e 2 ) 2 

Since the journal is precessing and not rotating, every point in 
the journal has a velocity equal to ern. The force in equation (2-38) 
therefore appears as a damping coefficient times a velocity acting 
in the direction opposite the journal motion. 

The equivalent bearing damping is: 

P _ yRL 3 T r (2- 4 °) 

0 “ — 7 

2c 3 (1 - e 2) 3/2 
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For the uncavitated film the limits of integration in 
equation (2-36) become: 


= 0, 0* = 2ir 

4U 


(2-41) 


Integrating and evaluating at those limits yields force 
components given by: 


F = 0 
r 


(2-42) 


F = -yRL 3 Treu (2-43) 

6 c 3 ( 1 - e 2 ) 3/2 

It is therefore evident that a complete fluid film does not 
produce an equivalent bearing stiffness but doubles the damping 
of the cavitated film. 

Although the equations for the bearing characteristics were 
derived for a plain bearing with no circumferential oil groove 
they are applicable to other bearing configurations. For ex- 
ample Figure (2-5 a) represents a plain bearing with circumferential 
oil groove and full end leakage. 

The total length of the bearing, L, corresponds to the 
length of the plain bearing with no oil groove. The bearing 
in Figure (2.5a) consists of two plain bearings without an oil 
groove whose length is L/2. Thus the bearing parameter L in 
the equations may be replaced by L/2 and the equations multiplied 
by 2 to obtain the total effect of the two half bearings. The 
net effect is to decrease the maximum pressure by a factor of 2: 
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(2-44) 



Similarly the damping and stiffness values are decreased by 
a factor of 4: 

1 (2-45) 

4 

The bearing represented in Figure (2.5b) is a plain bearing 
with circumferential oil groove and end seals to prevent end leak- 
age. If there is no end leakage the boundary conditions are; 




(2-46) 


and the net effect leaves the pressure and bearing characteristic 
equations unchanged. 



Figure 2^5 a. Axial Pressure Distribution of Bearing With 

Circumferential Oil Groove 
b. Axial Pressure Distribution of Bearing With End 
Seals and Circumferential Oil Groove 
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The bearing equations derived in this section are summarized 
in Table (2-1) , Also included in the table are the equations for 
pure radial squeeze motion. For this type of operation $ = o and 
results from a purely unidirectional load on the journal. The 
radial and tangential force components are derived from equation 

o 

(2-30) where only the term containing e in the integral is re- 
tained. The pressure equation is also modified to include only 

O 

the e term. The maximum pressure occurs at 0 = 7 r for all values 

of journal eccentricity. Examination of the pressure equation 

reveals that the hydrodynamic pressure is positive only in 
it 3tt 

the region i0' = — to These values of 0' are the limits of 

integration in equation (2-30) for the cavitated film. 

The table also shows that for. purely radial motion no bearing 
stiffness is obtained in either the cavitated or uncavitated 
bearing. Thus if this type of motion exists retainer springs 
must be included to provide support flexibility. 

For the case of circular journal precession, the table 
shows the stiffness and damping of the cavitated film and damp- 
ing of the uncavitated film remain essentially constant for low 
eccentricity ratios. As the eccentricity ratio increases above 
0.4 there is a rapid increase in these properties and they approach 
infinity as e approaches 1. This variation of stiffness in the 
cavitated film is very important. As the eccentricity becomes 
large the support becomes more rigid with a corresponding in- 
crease in the rotor critical speed. If the rotor critical speed 
is increased above the operating speed, the phase angle between the 
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Table 2-1. Summary of Equivalent Stiffness and Damping Coefficients 
for Squeeze Film Damper Bearings. 
















rotor unbalance vector and amplitude vector becomes less than 
90°. When this condition occurs the force transmitted through 
the support structure will always be greater than the unbalance 
load. With an uncavitated film this problem does not occur because 
no hearing stiffness is generated. To obtain the stiffness re- 
quired to stabilize a rotor (see Chapter 3) it is necessary to 
use retainer springs in the support bearings. 

One of the most significant parameters affecting damper per- 
formance is the length to clearance ratio. The stiffness and 
damping coefficients vary as (L/C)^ and therefore either doubling 
the bearing length or decreasing the clearance by 1/2 will in- 
crease, the coefficients by a factor of 8. 
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CHAPTER 3 


ROTOR-BEARING STABILITY AND -STEADY-STATE ANALYSIS 
3 . 1 ROTOR-BEARING STABILITY 

After Jeffcott's [5] analysis in 1919 of the single mass 
flexible rotor on rigid bearings, manufacturers began producing 
light, flexible rotors operating above the first critical speed. 
However , some manufacturers encountered severe operating diffi- 
culties with some of their designs. These machines underwent 
violent whirling while running above the critical speed and often 
failed . 

Experimental and analytical investigations by Newkirk and 
Kimball [6] [7J revealed that the whirl instability was not caused 
by unbalance in the rotor, but by internal shaft effects such as 
internal friction. Kimball theorized that forces normal to the 
plane of the deflected rotor could be produced by alternating 
stresses in the metal fibers of the shaft. In light of this theory, 
Newkirk concluded also that the same normal forces could be pro- 
duced by shrink fits on the rotor shaft. By incorporating these 
forces in Jeffcott’s model Newkirk showed that the rotor was un- 
stable above twice the rotor critical speed. 

Further investigation by Newkirk showed cases of rotor in- 
stability which were not produced by shaft effects but by effects 
in the journal bearings. [8] One cause of journal bearing in- 
stability was later shown to be due to lack of radial stiffness 
in the bearing and the instability occured at twice the rotor 
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critical speed. These instabilities were especially common in 
lightly loaded rotors and larger bearing loads tended to promote 
stability. The effect of the larger loads is to cause cavitation 
of the fluid film which results in a radial stiffness component 
of the bearing forces being produced. [9], [10], [11] 

In 1965 Alford reported on the effects of aerodynamic forces 
on rotors [12] . He showed that these forces couple the rotor 
equations of motion and can produce instability. He also noted 
that labyrinth seals and balance pistons also produce forces that 
can promote instability. 

Recent investigators including Gunter, Kirk and Choudhury 
[13] [14] [15] have analyzed the effects of support flexibility 
and damping on reducing rotor instability produced by the forces 
just described. As a result they have derived stability criteria 
for determining the necessary support characteristics. 

One of the most general methods for determining rotor 
stability is to derive the characteristic frequency equation of 
the system. The stability is given by the roots of this equation. 
The real part of the root corresponds to an exponentially in- 
creasing or decreasing function of time. Thus a positive real 
part indicates instability whereas a negative real part in- 
dicates a stable system. This type of stability analysis of a 
rotor-bearing system therefore requires that the characteristic 
equation be known. This equation is not always easy to obtain. 

The characteristic equation is derived from the homogeneous 
second order differential equations of motion of the system [15]. 
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By assuming solutions of the form 

x^ = A^e^ ^ i * 1, 2 , ...» n 

and differentiating, the equations are substituted back into the 
equations of motion. This produces a matrix known as the character- 
istic matrix. The determinant of this matrix gives the characteristic 
equation, a polynomial of degreee 2n in A, where n is the number 
of degrees of freedom of the system. 

The computer program SDSTB [16] was used to produce the 
stability maps shown in this chapter. The program calculates 
the characteristic equation for a three-mass symmetric flexible 
rotor mounted in journal bearings and supported in squeeze film 
damper bearings. The rotor-bearing model is shown in Figure (3-1). 

The rotor is assumed to remain stationary in the axial direction 
so the rotor has six degrees of freedom and the characteristic 
equation is therefore of degree twelve. The characteristic matrix 
is shown in Figure (3-2) . The determinant of this matrix gives 
the characteristic equation. The unknown variable in this 
equation is A, the natural frequency of the system. An examina- 
tion of the characteristic matrix shows that the coefficients 
of A are functions of the rotor and bearing properties as well as in- 
ternal shaft friction, absolute rotor damping and aerodynamic cross 
coupling.* The natural frequencies and stability of the system 
are found by finding the roots of this equation. 

The journal and support bearing characteristics can either be 
inserted directly as linear coefficients or they may be calculated 
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Figure 3--2 Characteristic Matrix For Three-Mass Model 

Including Aerodynamic Cross-Coupling Internal 
Shaft Friction and Absolute Shaft Damping 



in the program from the bearing parameters by solving for the 
equilibrium positions of the journal and support. These charac- 
teristics are non-linear functions of the journal eccentricity. 

The stability maps in this chapter were produced with the linear- 
ized journal and support bearing characteristics given as input 
data to the program. The assumption of linear bearing charac- 
teristics is usefulbbecause for low eccentricity the characteristics 
do not vary greatly with changes in eccentricity. This assumption 
allows a large savings in computer time. If the non-linear charac- 
teristics are calculated, the amount of computer time increases 
because an iterative procedure is used to find the equilibrium 
position. 

As an example of how a stability map is produced, consider 


the following system: 

ROTOR CHARACTERISTICS 
ROTOR WEIGHT 
JOURNAL WEIGHT 
SUPPORT WEIGHT 
SHAPT STIFFNESS 
SHAFT DAMPING 
INTERNAL DAMPING 
ROTOR SPEED 

BEARING CHARACTERISTICS 
k 

xx 

k 

yy 

c 

XX 


675 lbs 

312 lbs (each) 
15 lbs (each) 
280000 lb /in 

.10 lb-sec/in 
0.0 lb-sec/in 
10000 RPM 

1.287 x 10 6 lb/in 
1.428 x 10 6 lb/in 
1200 lb-sec/in 
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TJ 

yy 


1290 lb-sec/in 
0.0 lb/in 
0.0 lb-sec/in 


k = k 
xy yx 

C = 0 
xy yx 

Two values of aerodynamic cross coupling were selected, Q = 
20000 lb/in and Q » 100,000 lb/in. For each value of Q, several 
values of support stiffness were selected ranging from 50,000 
lb/in to 500,000 lb/in. For each value of support stiffness a 
range of support damping values from 0 to 10000 lb-sec/in was 
used. Using this method a stability contour was found for a 
given value of aerodynamic cross coupling and support stiffness. 

The rotor and bearing characteristics remained unchanged. 

Figures (3-3) and (3-4) show the stability maps for the 
above system for the two values of aerodynamic cross coupling. 

There is an intermediate range of support damping values for which 
the sytem is stable for a given value of the support. stiffness. 

As the stiffness is increased the system becomes less stable. 

With Q = 20000 lb/in the optimum amount of damping ranges from 
500 to 2500 lb-sec/in as the stiffness increases from 50000 to 
500000 lb/in. For damping less than 100 Ib-sec/in the system is 
unstable for all values of stiffness. The same is true if the 
damping exceeds 10000 lb-sec/in. 

For Q 100000 lb/in the optimum dampipg is 1000 lb-sec/in 
and does not shift over the stiffness range selected. When the 
stiffness reaches 250000 lb/in the system is unstable for all values 
of damping. 
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Figure 3-4 Stability of a Flexible Rotor With Aerodynamic Cross 
Coupling (Q = 100,000 lb/ln., N = 10,000 RPM) 


The stability calculation depends upon the accuracy of the 
root solving technique applied to the characteristic equation. 

For symmetric bearing and support characteristics repeated roots 
of the equation occur and many root solving routines break down 
under this condition. Although no cases have been encountered 
where the root solving routine has failed, it is possible that 
some combinations of rotor-bearing properties might cause this 
to happen. However because root solving routines are generally 
easy to obtain it would be easy to replace the one currently 
in SDSTB if such a situation arose. 

The amount of computer time required to solve the character- 
istic equation depends upon the root solving technique and the 
order of the characteristic equation. Many studies do not re- 
quire extensive stability maps and it is only necessary to de- 
termine whether the system is stable and not how stable. In 
these cases application of the Routh stability criteria [15] gives 
the required information without solving for the roots of the 
characteristic equation. This results in a savings of computer 
time. The option of using only the Routh stability criteria^ 
is available in SDSTB. 

The method of determining the stability of the system from 
the characteristic equation becomes less practical when the order 
of the system is large. The elements of the characteristic matrix 
must be found from the equations of motion and unless the deter- 
minant is found by a computer routine, the characteristic equation 
must be expanded by hand. Therefore complicated systems may re- 
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quire a prohibitive amount of formulation time. 

Recent research in transfer matrix and finite element 
techniques for determining the stability and natural frequencies 
of rotor bearing systems has been directed toward overcoming 
these difficulties [17] [18]. However, because iterative pro- 
cedures are required in the solution, higher order modes may be 
expensive to obtain from the standpoint of computer time. The 
type of method used will depend on the amount of computer funds 
available and the availability of programs using the various 
techniques . 

3.2 STEADY-STATE ANAJLYSIS 

The steady-state stability maps just discussed provide in- 
formation on the support characteristics needed to promote 
stability in a given rotor-bearing system. There remains the 
problem of relating these characteristics to the actual support 
bearing. The squeeze bearing equations derived in Chapter 2 in 
rotating coordinate's are used to determine the preliminary 
bearing design. As noted in Chapter 2, these equations were de- 
rived assuming steady-state circular, synchronous precession of 
the journal. 

The bearing characteristics, stiffness, damping and pressure 
are functions of the amplitude of the journal orbit, fluid vis- 
cosity and bearing geometry. The addition of oil supply grooves, 
end seals and cavitation affect the bearing characteristics. 

The steady-state equations have been programmed on a digital 
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computer. This program, SQFDAMP, analyzes three basic bearing 
configurations : 

1. Plain bearing, no oil supply groove or end seals. 

2. Bearing with oil supply groove but without end seals. 

3. Bearing with both oil supply groove and end seals. 

Both cavitated and uncavitated fluid films can be analyzed. A 

listing of the program with a description of the input data re- 
quirements and sample output are contained in Appendix A. 

The program calculates the bearing characteristics and plots 
them as functions of the journal eccentricity ratio, e. By 
varying the bearing parameters the designer is able to determine 
the bearing characteristics and select a bearing configuration 
that will provide the stability requirements of the system under 
consideration. 

Figures (3-5) - (3-7) show the characteristics for a bearing 
being considered for the 675 lb rotor system described earlier. 

The bearing has an oil supply groove and end seals, and the fluid 
film is assumed to be cavitated. The bearing parameters are, 
length, 1.0 inches, radius, 1.2 inches and fluid viscosity 10 
microreyns. For the case where Q ■ 20000 lb/in it was determined 
that the optimum support damping is about 500 lb-sec/in and 
the support stiffness should be less than 100000 lb/in. Because 
it is desirable to keep the eccentricity ratio of the journal 
low, Figures (3-5) and (3-6) reveal that this bearing will provide 
the necessary stiffness and damping characteristics with a clear- 
ance of about 4 mils at an eccentricity ratio of e = .10 to *20. 
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Figure 3-5. Damping Coefficient for Squeeze Film Bearing 

With Cavitated Film - End Seals and Oil Supply 
Groove Included. 
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STIFFNESS, K0 - (LB/IN 1 


SQUEEZE FILM DAMPER 
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§ CAVITATION 
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Figure 3-6. Stiffness Coefficient for Squeeze Film Bearing 
With Cavitated Film - End Seals and Oil Supply 
Groove Included. 
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Figure 3-7. Maximum Pressure for Squeeze Film Bearing with 
Cavitated Film - End Seals and Oil Supply 
Groove Included. 
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This corresponds to a journal orbit of 0.4 to 0.8 mils amplitude. 
The maximum hydrodynamic pressure in the bearing is about 100 
psi for this clearance. If the fluid film cavitates, the re- 
sulting characteristics are shown in Figures (3-8) and (3-9) . 

A slightly larger clearance, 5,0 mils, will produce the optimum 
damping. However, because the uncavitated film does not produce 
an equivalent stiffness, retainer springs must be incorporated 
in the bearing. If the end seals are flexible the required 
spring rate may be obtained from them. 

One advantage of the uncavitated film is that if the journal 
eccentricity ratio should become very large, there is no rise in 
stiffness that could cause the system to become unstable or 
raise the critical speed above the operating speed. Figures (3-5) 
and (3-8) indicate that even at eccentricity ratios of 0.9 the 
damping value still remains acceptable. For the cavitated film 
at e - .9, the stiffness exceeds 2,000,000 lb/in and the system 
would be bordering on instability. Fox both films the maximum 
pressure exceeds 70000 psi. at e *= 0.9 and this large a rotating 
pressure field could result in bearing failure. 

If Q = 100000 lb/in, the bearing characteristic graphs 
reveal that for a cavitated film the clearance must be as small 
as possible because of the limitations on stiffness shown in 
Figure (3-4). With a 3.0 mil clearance, the stiffness is 250000 
lb/in at e = 0.23 and this stiffness will produce system in- 
stability. For Q = 100000 lb/in it is desirable to study other 
bearing lengths and radii to obtain a cavitated bearihg which pro- 
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Figure 3-8. Damping Coefficient for Squeeze Film Bearing 
with Uncavitated Film »■ End Seals and Oil 
Supply Groove Included 
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Figure 3-9. Maximum Pressure for Squeeze Film Bearing 
with Uncavitated Film - End Seals and Oil 
Supply Groove Included. 
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duces the required damping at low eccentricity ratios and also 
produces a stiffness of 100000 lb/in or less. For an uncavitated 
film a 3.0 mil clearance provides adequate damping. 

The steady-state bearing characteristic equations used in 
conjunction with the stability analysis based on steady-state 
motion provide an excellent means of determining bearing con- 
figurations. Using the methods just described, good preliminary 
designs can be obtained which can be more thoroughly analyzed. 

The steady-state analysis described here and the transient analysis 
described in the next chapter provide bearing design criteria 
which will eliminate the experimental testing of unsuitable de- 
signs. The costs of the analysis more than offset the experimental 
losses when designs are tested that result in bearing or machine 
failure. 
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CHAPTER 4 


TRANSIENT ANALYSIS 

4.1 INTRODUCTION 

There are a number of operating conditions in which the 
squeeze bearing journal does not orbit the bearing center in 
circular synchronous precession. These conditions can occur when 
there is a unidirectional load on the rotor or when there is a 
suddenly applied load such as the application of unbalance when 
blade loss occurs. Intermittent or cyclic forces transmitted 
to the machine from nearby equipment can also result in non- 
linear orbiting. Under these conditions the bearing stiffness 
and damping coefficients developed in Chapter 2 are no longer 
applicable and a time-transient analysis of the bearing is 
necessary to determine the squeeze bearing support's ability to 
restablize the system. 

The equation for the squeeze film damper fluid film forces, 
in fixed coordinates, equation (2-19 \ provides a useful means of 
determining the time dependent transient behavior of a rotor 
bearing system. The force equations have been programmed on a 
digital computer and combined with the journal equations of motion. 
The resulting program, BRGTRAN [19] tracks the journal motion 
forward in time under the influence of the bearing forces, journal 
weight and unbalance. The program is capable of including the 
effects of retainer springs and cavitation. 

Because both the bearing pressure equation and the journal 
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equations of motion must be integrated, the accuracy of the sim- 
ulation depends upon the numerical integration method used. Two 
integration methods, Adam-Bashforth-Moulton Predictor-Corrector, 
and 4th Order Runge- Kutta, are provided as options in BRGTRAN 
Although the 4th Order Runge-Kutta method is highly accurate, 
four functional evaluations of the pressure equation are required 
for each step in time. Even though the required time step size 
may be larger using the Runge-Kutta method, the total computer 

time required is still greater than other methods because of the 
many functional evaluations. 

The Adams-Bashf orth-Moulton method has been found to be 
sufficiently accurate for most cases run if the time step has 
been made sufficiently small, or about 0.01 cycles. At low 
eccentricities no problems are encountered in the integration 
process and the solutions are accurate for both methods. The 
changes in the bearing forces are relatively small from one time 
step to the next and even a very simple integration method such 
as the Modified Euler Method provides reasonable accuracy. However 
at high eccentricities the bearing forces change drastically with 
even a very small change in eccentricity. Even the more sophisti- 
cated integration schemes lose accuracy when the eccentricity is 
high unless the time step is made very small. The amount of 
computer time and core storage required for a small time step 
becomes prohibitive. This is especially true in light of the 
fact that at high eccentricities the validity of the short bear- 
ing approximation used in reducing the Reynolds equation is 
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doubtful. In using the short bearing approximation it was 
assumed that the pressure gradient in the tangential direction 
is small, and this assumption may be violated at high eccent- 
ricities. 

It was also shown in Chapters 2 and 3 that at high eccent- 
ricities the bearing stiffness becomes very large. This can 
result in system instability as indicated in the stability maps, 
Figures (3-3) and (3-4), or in raising the system critical speed 
above the operating speed causing large forces to be transmitted 
to the machine structure. For these reasons the design criteria 
of e<0.4 was established and it is unnecessary to use excessive 
computer time to obtain greater accuracy at higher eccentricites , 
The Information provided at these eccentricites is very useful 
in showing trends in the ability of a bearing to perform adequate- 
ly and should be used with this restriction kept in mind. Also 
recent transient analysis has been performed using hybrid com- 
puter simulation thereby avoiding the difficulties inherent in 
numerical integration. [20] 

4.2 ANALYSIS 

Using a time transient bearing program, a design engineer can 
make an analysis of the bearing effects without resorting to 
either a complete time-transient analysis of the entire rotor- 
bearing system or a costly experimental program during the 
preliminary design stage. During later design stages when the 
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bearing ( configuration has been tentatively fixed, a more com- 
plete theoretical and experimental analysis of the entire 
system may be performed. The bearing force calculations have 
also been incorporated as a subroutine in a computer program 
which analyzes the transient behavior of certain rotor-bearing 
models [21]. By using these programs to determine the bearing 
parameters experimental verification of the bearing effects 
can be performed with more assurance that the design is feasible, 
and costly and time consuming machine prototype failures can be 
reduced. 

The computer program BRGTRAN was recently used as part of 
an analysis of an existing turbomachine which had suffered 
frequent bearing failures. The manufacturer had decided to use 
a squeeze film damper bearing to reduce the vibration amplitudes 
at the failing bearings. Without performing a complete analysis 
of the bearing effects, an experimental program was initiated 
where various bearing configurations were installed on a test 
machine. The damper bearings used did not dampen j out the vibrations 
and much time and money was lost during the project. 

The following discussion of the computer simulation of this 
system shows the effect of varying the rotor-bearing parameters 
including bearing length, unbalance, cavitation and retainer springs. 

In this study of the problem, the analysis made using BRGTRAN 
showed why the bearings used were unable to reduce vibration am- 
plitudes. Figure (4-1) shows the journal orbit in a cavitated 
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Figure 4-1. Unbalanced Rotor in Cavitated Squeeze Film 

Bearing - L = 0.45 IN. - Unbalance Eccentricity 
* 0.002 IN. - No Retainer Spring. 
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film during the first 10 cycles of transient motion. The bear- 
ing parameters are: 

LENGTH - 0.45 inches 

RADIUS - 2.55 inches 

CLEARANCE - 0.004 inches 

FLUID VISCOSITY - 0.38 X 10" 6 

JOURNAL WEIGHT --74 lbs. 

UNBALANCE ECCENTRICITY - 0.002 inches 
RETAINER SPRING RATE 



Kxx 

- 0 lb/in 


Kyy 

- 0 lb/in 

N 


- 16800 RPM 


In the transient orbit figures a standard right-hand co- 
ordinate system has been adopted with positive journal rotation 
in the counterclockwise direction. The .asterisk on the orbit 
represents the point where the maximum force is generated. 

The small dots represent timing marks denoting one revolution 
of shaft motion. These marks can be used to determine the re- 
lative location of the unbalance with respect to the amplitude 
in the x-direction. The timing mark sensor is assumed to be 
located on the positive x-axis, and the phase angle is measured 
from the x-axis to the timing mark in the clockwise direction. 

Returning to Figure (4-1) , it is seen that the journal very 

quickly spirals out to an eccentricity ratio e = 0.95, where a 

limit cycle is formed due to the non-linearity of the bearing 
forces. The maximum force transmitted through the support 
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structure is .7405 lbs, which is over 6 times the rotating un- 
balance load of 1181 lbs, as shown by the parameter TED 
Figure (4-1) - This large rotating force will eventually lead 
to bearing failure and is therefore undesirable. 

The phase angle between the maximum amplitude in the x 
direction and the timing mark is approximately 30° , and this 
indicates that the precession rate is less than the natural 
frequency of the bearing. This has been caused by the large 
stiffness developed in the bearing. Figure (4-2) shows that 
the stiffness is approximately 487,000 lb/in. The damping is 
given in Figure (4-3) as 71,5 lb-sec/in. 

One possible design change considered was to increase 
the bearing length. Figure (4-4) shows the effect of increas- 
ing the bearing length of 0,90 inches, all other parameters 
remaining the same. The journal still spirals outward, how- 
ever the limit cycle is produced at an eccentricity ratio, e * 
0.88, This results in a reduction in the force transmitted to 
the support from 7405 to 3371 lbs,, but it is still greater 
than the rotating unbalance load, and is undesirable since bear- 
ing failure will result. The phase angle has shifted from 
30° to 60° and this would have been accompanied by an increase 
in the force transmitted except that the stiffness and damping 
values have also changed with the net effect being a reduction in 
the transmissability. The stiffness and damping coefficients are 
shown in Figures (4-5) and (4-6) to be 880000 Ib/in and 22 lb- 
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Figure 4-3. Damping Coefficient for Squeeze Film Bearing 
of Figure (4-1). 
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Figure 4-5. Stiffness Coefficient for Squeeze Film Bearing 
of Figure (4-4). 
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Figure 4-6. Damping Coefficient for Squeeze Film Bearing of 
Figure (4-4) . 
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sec/in respectively. 

Figure (4-7) shows the first 10 cycles of motion for the 
bearing configuration of Figure (4-1) with the addition of re- 
tainer springs with a stiffness of 123000 lb/in. Only a very 
slight improvement in force transmission results, and the bearing 
is still unacceptable. 

The effect of changing the rotating unbalance load is 
shown in Figure (4-8) . The bearing configuration is the same 
as Figure (4-4) except the unbalance has been reduced by one 
half. The journal orbit has been greatly reduced and the force 
transmitted to the support structure is reduced to only 80% of 
the unbalance load. The journal is now precessing about a point 
offset from the bearing center. 

By adding :retainer springs to the bearing of Figure (4-8) 

the motion in Figure (4-9) results. The retainer spring rate 

is 123,000 lb/in. The journal is orbiting about the center of 

the bearing at an eccentricity ratio of e = 0.35. The journal 

motion is stabilizing more quickly than without the retainer 
springs and the journal is precessing synchronously as indicated 

by the small dots on the or bit which represents one cycle of 

motion. The transmitted force has been further reduced to only 

65% of the unbalance load. From the standpoint of producing 

a small amplitude orbit and attenuating the unbalance load such 

a bearing configuration is desirable. 

If the unbalance eccentricity is ag£in increased to 0.002 
inches, the bearing configuration of Figure (4-9) is no longer 


54 



SQUEEZE FILM BEARING- 
CAVITATED FILM 





HORIZONTAL 

esse 

sasro 

W 

r 

73.7 LBS 

N = 

16000 

RPM 

L 

— 

.450 IN 

R = 

2.550 

IN 

C 

r 

4.00 MILS 

MU = 

.302 

MICROREYNS 

PS 

.s 

0.00 PSI 

FMAX= 

7331.3 

LBS 

WX 

r 

0.00 L3S 

WY = 

0.00 

LBS 

FU 

s 

1181.31 LBS 

EMU - 

.50 


KRX 

n 

123000 LB/IN 

KRY = 

123000 

LB/IN 

TRO 

— 

6.20 

PMAX= 

12510.22 

PSI 



Figure 4-7. Unbalanced Rotor in Cavitated Squeeze Film Bear- 
ing - L = 0.45 IN. - Unbalance Eccentricity « 
0.002 IN. - Retainer Spring Stiffness, KR = 
123,000 LB/IN. 
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Figure 4-8. Unbalanced Rotor in Cavitated Squeeze Film Bear- 
ing - L = 0.90 IN. - Unbalance Eccentricity = 
0.001 IN. - No Retainer Spring 
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Figure 4-9. Unbalanced Rotor in Cavitated Squeeze Film Bear- 
ing - L *= 0.90 IN - Unbalance Eccentricity = 
0.001 IN. - Retainer Spring Stiffness, KR = 
123,000 LB/IN. 
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acceptable. The resulting motion is shown in Figure (4-10) 
and a large amplitude limit cycle with large force transmission 
again occurs. However if the fluid film does not cavitate the 
bearing performance improves and is marginally acceptable as 
shown in Figure (4-11). The journal is orbiting at e = 0.65 
and the transmitted force is only 5% less than the un- 
balance load. 

Although this particular analysis includes only the journal 
weight and the unbalance loading, it shows the usefulness of 
this program in determining the bearing effects with different 
bearing parameters. The ability to perform the analysis without 
extensive preliminary experimental work provides a great savings 
in time and money. By systematically varying the bearing para- 
meters design guidelines are established. 

For instance Figures (4-12 - (4-15) shows the effect of 
varying the unbalance on a 675 lb. journal operating in a squeeze 
bearing with a 7 mil clearance. The unbhlance eccentricity is 
increased from 1.75 to 3.5 mils with an accompanying increase 
in the journal amplitude and the force transmitted to the support 
structure. Although the exact unbalance may not be known pre- 
cisely for a given rotor, a design estimate can be made based on 
the effect of a suddenly applied known unbalance due to blade loss 
or loss of chemical deposits from the blade surfaces. Prior to 
the sudden unbalance the rotor is assumed to be perfectly balanced. 
The ability of the bearing to reduce the amplitude to tolerable 


58 



SQUEEZE FILM BEARING 
CAVITATED FILM 
HORIZONTAL 


N 

r 

73.7 LB5 

L 


.900 IN 

C 

r 

4.00 MILS 

PS 

- 

0.00 PSI 

NX 

r 

p.00 LBS 

FU 

— 

1181.91 LBS 

KRX 

r 

123000 LB/IN 

TRO 

— 

3.08 


N = 16800 RPM 

R = 2.550 IN 

MU = .382 MICR0REYN5 

FMAX= 3641.9 LBS 
NY = 0.00 LBS 

EMU = .50 

K.RY = 123000 LB/IN 

PMAX= 3334.52 PSI 



Figure 4-10. Unbalanced Rotor in Cavitated Squeeze Film Bear- 
ing - L « 0.90 IN. - Unbalance Eccentricity - 
0.002 IN. - Retainer Spring Stiffness, KR * 
123,000 LB/IN. 
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Figure 4-11. Unbalanced Rotor in Uncavitated Squeeze Film Bear 
ing - L * 0.90 IN. - Unbalance Eccentricity = 

0.002 IN. - Retainer Spring Stiffness, KR = 

123,000 LB/IN. 
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Figure 4-12. Vertical Unbalanced Rotor in Squeeze Film Bear- 
ing - Effect of Unbalance Magnitude - Unbalance 
Eedentricity .** 1.75 Mils 
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Figure 4-13. Vertical Unbalanced Rotor in Squeeze Film Bear- 
ing - Effect of Unbalance Magnitude - Unbalance 
Eccentricity * 2.10 Mils. 
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Figure 4-14. Vertical Unbalanced Rotor in Squeeze Film Bear- 
ing - Effect of Unbalance Magnitude — Unbalance 
Eccentricity = 2.45 Mils. 
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Figure 4-15. Vertical Unbalanced Rotor in Squeeze Film Bear- 
ing - Effect of Unbalance Magnitude - Unbalance 
Eccentricity = 3.50 Mils. 
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levels until the machine is shut down or its operating conditions 
changed can be determined. Note the shift in phase angle from 
180° to 90° as the eccentricity increases. 

It has been shown that retainer springs help center the 
journal and reduce the vibration amplitude. They may also be 
used to prevent oil leakage from the end of the bearing if they 
are of the 0 ring type. As noted in Chapter 2, the uncavitated 
film provides no equivalent stiffness when the journal is operat- 
ing in synchronous precession about the bearing center. In this 
case the use of a retainer spring to provide a restoring force 
in the bearing is necessary. Although an increase in stiffness 
results in centering the journal in the bearing, (see Figures (4-16) 
(4-19) the magnitude of the transmitted force is minimized for some 
intermediate value of stiffness which depends on the bearing para- 
meters and loading. Also the ability of the journal to quickly 
return to a stable, steady state operating condition -is impaired 
with increasing stiffness. Both of these operating conditions 
were indicated by the stability maps in Chapter 3. 

After the preliminary bearing design, a more thorough 
analytic study may be made by incorporating the damper bearing 
effects into a program which includes the dynamic effects of the 
entire system. Although such programs are not readily available 
for many complex systems, one has been developed for a three mass 
rotor in journal bearings on a squeeze film bearing support as 
noted earlier. The rotor bearing model is shown in Figure (3-1)* 
This program may be used to calculate the ability of squeeze 
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Figure 4-16. Horizontal Unbalanced Rotor in Squeeze Film Bear- 
ing - Effect of Retainer Springs - Retainer 
Spring Stiffness, KR = 0 LB/ IN. 


66 




SQUEEZE FILM BEARING- 
CAVITATED FILM 






HORIZONTAL 

case vo. 

3ZS7U13 

w 


G75.0 

LBS 

N = 

10500 

RPM 

L 

— 

2.000 

IN 

R = 

3.500 

IN 

C 

Z2 

15.00 

MILS 

MU = 

2.490 

MICROREYNS 

PS 

ZZ 

0.00 

PSI 

FMAX= 

1054.1 

LBS 

wx 

— 

0.00 

LBS 

WY = 

0.00 

LBS 

FU 

— 

351.40 

LBS 

EMU = 

.03 


KRX 

— 

50000 

LB/IN 

KRY = 

50000 

LB/IN 

TRD 

— 

1.11 


PMAX= 

39.41 

PSI 



Figure 4-17. Horizontal Unbalanced Rotor in Squeeze Film Bear- 
ing - Effect of Retainer Springs - Retainer 
Spring Stiffness, KR = 50,000 LB/IN. 


67 




SQUEEZE FILM BEARING 
CAVITHTED FILM 






HORIZONTAL 

Oftt MS. 

s&rm 

w 

— 

675.0 

LBS 

N = 

10500 

RPM 

L 

- 

2.000 

IN 

R = 

3.500 

IN 

C 

« 

15.00 

MILS 

MU = 

2.430 

MICR0REYNS 

PS 

- 

0.00 

PSI 

FMAX= 

336.7 

LBS 

wx 

— 

0.00 

LBS 

WY = 

0.00 

LBS 

FU 


351.40 

LBS 

EMU = 

.03 


KRX 


100000 

LB/IN 

KRY = 

100000 

LB/IN 

TRD 

r 

.38 


PMAX- 

6.10 

PSI 



Figure 4-18. Horizontal Unbalanced Rotor in Squeeze Film Bear- 
ing - Effect of Retainer Springs - Retainer 
Spring Stiffness, KR = 100,000 LB/IN. 
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film bearings to stabilize a multi-mass rotor. In adidtion, in- 
formation may be obtained to verify the stability maps obtained 
by other means (see Chapter 3). Because the bearing is initially 
designed using the criteria derived from a stability analysis 
such verification is useful in judging the overall worth and 
limitations of the analytical design process. Using these analy- 
tical methods leads to a more efficient testing program because 
unacceptable bearing designs are eliminated before testing begins. 
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CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 

The problem of rotor stability is of current importance 

because of the high speeds and complex dynamics of modern rotor 

bearing systems. Damped flexible supports have a great effect 

on the ability of a system to suppress unstable whirl. There- 
fore there is a need for methods of predicting rotor instability 

and obtaining bearing designs that provide the necessary support 

characteristics . 

5.1 PREDICTING ROTOR INSTABILITY 

Methods for determining the stability of rotor bearing systems 
include: 

1. Using the rotor^bearing system equations of motion in- 
cluding the effects of aerodynamic forces, shaft damp- 
ing and internal friction to obtain the system 
characteristic equation. The roots of this equation 
show the stability and natural frequencies of the 
system. 

2. Using finite element and transfer matrix methods to 
obtain the system stability and natural frequencies. 

3. Applying Routh stability criterion to the system 
characteristic equation. 

The third method does not yield information on relative 
stability but only determines whether or not the system is absolutely 


71 



stable. The first two methods provide relative stability in- 
formation that can be used to plot stability contours as funct- 
ions of the support characteristics. The stability maps pre- 
sented show that for a given value of support stiffness there 
is a range of damping values which will stabilize the system. 

If the support stiffness becomes too large, the system will be 
unstable for all values of support damping. The optimum 
stiffness and damping values for a particular system depend 
upon the rotor-bearing properties and the nature and magnitude 
of the forces acting on the system that produce instability. 

5.2 DETERMINING THE STIFFNESS AND DAMPING COEFFICIENTS OF 
THE SQUEEZE FILM DAMPER BEARING 

The assumption of steady state circular synchronous pre- 
cession of the journal and the use of a rotating coordinate 
system allow the bearing forces to be equated to equivalent 
stiffness and damping forces. This establishes stiffness 
and damping coefficients for the bearing. These coefficients 
are functions of the bearing geometry, the use of end seals and 
cavitation of the fluid film. The coefficients obtained from 
the steady state bearing analysis are compared with the values 
from the stability maps to determine a bearing configuration 
that promotes system stability. 

The stiffness and damping coefficients are non-linear 
functions of the journal eccentricity. However for e<0.4 the co- 
efficients do not change appreciably with changes in eccentricity. 
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For values of e>0. 4 a rapid increase in the coefficients 
occurs. The high stiffness developed can cause system in- 
stability or raise the system critical speed above the operat- 
ing speed resulting in force transmissabilities greater than 
1. For these reasons a design criteria of e<0.4 has been es- 
tablished. If the fluid film does not cavitate a radial stiff- 
ness is not developed and must be supplied by retainer springs. 

5.3 TRANSIENT ANALYSIS 

The bearings designed using steady state analysis are 
further analyzed using transient response programs. The motion 
of a system under the influence of unbalance and other external 
forces is monitored. Effect of retainer springs to preload the 
bearing can be determined and the bearing design further refined. 

The accuracy of transient response programs are dependent 
on the accuracy of the numerical integration methods employed, 

At highueccentricities very small integration step sizes are re- 
quired to retain high accuracy in the solution because of the 
rapid variation in the bearing forces, The cost of obtaining high 
accuracy at high eccentricity does not justify using small time 
steps or complicated integration methods requiring many functional 
evaluations because optimum bearing design requires operation at 
low eccentricities. The information provided by simpler inte- 
gration methods is sufficient to indicate trends in bearing operation 
at high eccentricities. 
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5.4 ADVANTAGES OF BEARING SIMULATION 


The analytic simulation of the squeeze film damper bearing 
eliminates many bearing designs that would result in bearing 
or machine failure in a test installation. Because the cost of 
constructing and instrumenting a test rig is very high, prevent- 
ing the failure of these machines is important. The time involved 
in manufacturing and testing bearings is also great and the elim- 
ination of unsuitable designs by analytic procedures results in 
a substantial savings in time and money. 

A good test program is essential, however, to determine the 
actual rotor-bearing response under various conditions. The ana- 
lytic simulation provides a means of interpreting the test data. 

Often the nature of the actual system excitation is unknown and 
the actual system response must be used to infer the nature of these ex- 
citations. Where it is possible to systematically vary the ex- 
citation the experimental results provide 3 . check on the accuracy 
and limitations of the analytic simulation. 

5.5 LIMITATIONS OF ANALYTICAL INVESTIGATIONS 

In any analytical investigation it is very important to know 
the assumptions made in analyzing the problem. In deriving the 
bearing equations for this study several assumptions were made. 

To obtain the Reynolds equation from the Navier-Stokes equations 
it was assumed that: 

1. The viscosity is constant 

2. The flow is steady state 

3. The fluid inertia terms are negligibly small 
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4. The density is constant 

5. There is no flow in the radial direction 

6. There is no pressure gradient in the radial direction 
The Reynolds equation obtained was modified by using the 

short bearing approximation. The assumption made was that the 
pressure gradient in the tangential direction is small and when 
multiplied by h^ it is very small in comparison to other terms 
containing only h. This assumption is valid only if the tan- 
gential pressure gradient is small, and at high eccentricity 
ratios this assumption may not be valid. For this reason the 
design criterion is that the eccentricity ratio be less than 0.4. 

In addition the short bearing approximation is valid for < 0.5. 

For ratios greater than this the axial pressure gradient is not 
large compared to the tangential gradient. For > 0.5 either 
the long bearing approximation or finite bearing techniques 
should be used. 

The conditions under which cavitation occurs must be modified 
in light of the current experimental results. In this study cav- 
itation conditions wete assumed to be the same in squeeze bearings 
as in journal bearings. 

5.6 RECOMMENDATIONS FOR FUTURE RESEARCH 

1. Construct an experimental rotor with squeeze damper 
supports to provide data to verify the analytic squeeze 
bearing model, 

2. Conduct analytic and experimental research into the con- 
ditions under which the squeeze bearing fluid film cavitates 
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and to determine how cavitation propagates through 
the film. 

3. Investigate the heat transfer characteristics of the 
bearings and provide modifications to the bearing pro- 
grams to account for variable viscosity of the lubri- 
cant. 
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APPENDIX A 


DESCRIPTION OF PROGRAM SQFDAMP 

This program analyzes the stiffness, damping and pressure 
characteristics of the squeeze film damper bearing. Three bearing 
configurations may be analyzed: 

1. Plain bearing without end seals or circumferential oil supply 
groove. 

2. Bearing without end seals but with circumferential oil supply 
groove. 

3. Bearing with both end seals and circumferential oil supply 
groove. 

In addition the fluid film may be assumed to be cavitated or un- 
cavitated. If cavitated the film extends from 0 = y to 0= If un- 

cavitated it extends from 0 = o to 0 = 2ir. The 0 is measured from the 
line of centers of the bearing in the direction of journal precession. 

The evaluation of the bearing characteristics assumes that the journal 
precesses synchronously about the bearing center in a circular orbit. 

The following is a description of the program input data: 

CARD 1 -80 column free-field comment card. 

CARD 2 -80 column free-field comment card. 

CARD 3 -Namelist/BRGTYPEi/ TYP, CAV, PS 

TYP - 0 for bearing type 1 (see above) 

- 1 for bearing type 2 (see above) 

- 2 for bearing type 3 (see above) 

CAV - 0 for uncavitated film 

- 1 for cavitated film 
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CARD 4 


CARD 5 
CARD 6 
CARD 7 


Sample input 
CARD 
1 
2 

3 

4 

5 


PS - oil supply pressure, psi. 

Nameli s t / BEARING/ L,R, MU, N 
L - Bearing length, in. 

R - Bearing radius, in. 

MU - Lubricant viscosity, microreyns 
N - Rotor speed (journal precession rate), RPM 
Namelist/ECRATIO/ES ,EF 

ES - initial journal eccentricity ratio ES>0 
EF - final journal eccentricity ratio EF<1 
Namelist /CLEARNC/ C(I), NC 
G(I) - clearance 

NC - Number of clearance values 
Namelist/PLOTSEM/ CS, PC, PK, PP 

CS - Plot control, .T. if plot desired, otherwise .E. 
PC - .T. if damping plot desired, otherwise .F. 

PK - .T. if stiffness plot desired, otherwise .F. 

(if CAV=0, PK = .F.) 

PP - .T. if pressure plot desired, otherwise .F. 

data. 

DATA 

SAMPLE DATA FOR SQFDAMP 
1 MARCH 1973 

$BRGTYP TYP = 0 , CAV = 1, PS = 0.0$ 

$BEARING L = 0.90, R = 2.55, MU = 0.382, 

N - 16800.0$ 

$ECRATI0 ES = 0.1, EF = 0.9$ 
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6 $CLEARANC C(l) = .003, C(2) = .004, C(3) # .005 
C(4) - .006, NC - 4$ 

7 $PLOTSEM CS = .T., PC = .T., PK = .T., ?P = .T.$ 
The following is a listing of the program and a sample output. 
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P_R_Q.GRAM__S.QEO AMP1I N.F.U.I ,.0_U \ P UJ, T_A_PE5 = IliP_yTjJ.APE.6jL0yiP.Ulj UPJU1 

C 

_c VER_SI0N_- 13 MAR CH 19 73 

C 

R£A L Uj.MU , K 0 , N __ 

LOGICAL CRING,CS,PC,PK,PF 

INT EGER_C„0!4 E HT 1 .C0HENT2 .T Y P 

INTEGER CAV 

OIMENSI.fi> £0(135) , C 0 ( 1 3 5 ) , K.0 ( 1 3_5 ) j P M A X_(1.3 5. ) , T H El A M_ ( 1 3_5)_ ,C (.5 ). j 

1C0MENT1 (8) , CCMENT2 (8) 

N A .HE LJ S.T./ BRGTYPE/ T_YP»CAV, FS 

NAMELIST/BEARING/ L,R,MU,N 

N.A M E.L I SJ /E.Q RATIO /_E S jEE 

NAMELIST/CLEARNC/ C,NC 

.NA M EL IS I/.PL.O IS E M/_ C_ S , PC. , P K j £P . . 

COMMON NEO 

Q : 

C READ DATA 

JL : 

K=1 

JaJLtt. RE AD (5 f 1 ) C0MENT1 

1 FORMAT ( 8 A10 ) 

I E.( E.CEj.5.)_9 999 , 2 . ; 

2 READ (5 » 1 ) C0MENT2 

_R£AD.(.5j.6RG.T YPE) 

READ (5 , 6EARING) 

RE AD.(5j ECRATIO) . , 

READ (5»CLEARNC) 

READ (5.j pL 0T.S.E Ml . . 

IF (.NOT. CS) GOTO 3 

JLEJLK-tEQ »..0.)__GO.IQ_ 3 

CALL CALCOMP (1) 

-C „ — — ~ — — — — — 

C SET UP NEW BLOCK ANG MOVE ORIGIN 

__C ALL__PL Qj 0_*_ - 

C 

JL 

C 

_J3 WRITE <6 rkl _ - 

4 FOR HAT ( 1H1, 35X, 27H****SCUEEZE FILM DAMFER*** 4 /) 

_WRI.JE <6 ,5_> _ 

5 FORMAT (IX, 

i5B.HT.HIS_ PROGRAM. ANALYZES ..THE. STIFFNESS DAMPING AND. PRESSURE „/JUb_ 

258H CHARACTER 1ST ICS CF THE SQUEEZE FILM CAMPER BEARING, THREE /IX, 

_ 3 5 8 H E E AR I MG — C C NF I G U R AT ION S MAY BE ANALYZED- _____ /XX$ 

458H 0 - PLAIN BEARING MTKCUT END LEAKAGE SEALS CR CIRCUM- /IX, 

... 558H FERENTIAL CIL SUPPLY GROOVE /lX f 

658 H 1 - BEARING WITH CUT END LEAKAGE SEALS BUT WITH CIRCUM- /IX, 

_ Z58H FERENTIAL CIL _ SUPPLY. GROOVE . _ ____ . 

858H 2 - BEARING WITH BOTH END LEAKAGE SEALS ANO SIRCUMFER- /IX, 

_ 958H ENTIAL OIL SUPPLY GROOVE _ /IX, 

158HIN Addition, the film may be assumed to be either cavi- /ix, 

___ 258HTATEC OR UNCAVITATED. IF CAV IT A TED THE FILM IS ASSUMED _ /IX, 
358HT0 EXTEND FROM THETA = FI/2 JO THET A-3FI/2 , WHERE THETA /IX, 

_ . ...458HIS MEASURED FROM THE LINE CF_ CENTERS IN THE DIRECTION OF _./iXj_. 

558HJ0URNAL PRECESSION. THE EVALUATION CF THE BEARING CHARAC- /IX, 
.658 KT ERISTICS ASSUMES. THAT. THE .JOURNAL FRECESSES SYNCHRONOUS- . /IX , 


SGFO 010Q__ 
SQF00110 
SQF 00 120 
SQF 00130 
SQF DO 140 
SQFD0150 
SQF DO 16 0_ 
SQFQ017 0 
_SQFD0180 
SQFD0190 
SQFOQ200_ 
SQF DO 210 
SCFD02ZO 
SQF DO 230 
SQFD024Q 
SQFD0250 
SGFDQ260 
SQF DO 270 
.SQF DO 28 0_ 
SQF DO 2 90 
_SQF0030 0 _ 
S OF 003 10 
_SQF0Q320 
SQF DO 330 
SQF0Q340_ 
SQF DO 350 
_SQFD0360 
SQF 00 37 Q 
SQF 00380 
SQF D 03 90 
SQFD040 0_ 
SQF 00410 
SQF 00420 
SGF00430 
SQF00440 
SQFD0450 
SQf D0460_ 
SQF00470 
„SCFD0480 ‘ 
SQF00490 
SQF 005 CO 
SQFOC510 
SQF DOS 20 
SQFD053Q 
SQF0C54Q. 
SCFD0550 
SQF DO 560 
SQFOC570 
SQFD0580 
SQFD0590 
SGF00600 
SQFDG610 
SQF00620 
SQFD0630 
SQF 00640 
SQFDu653 
SQF00660_ 
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__„ 6 - 


758HLY AECUT THE BEARING CENTER. /IX, 

J658HIHE FCLL OWING .IS „A._.D ESC ft IP.TI.C N. CFJHE, I NP.UX— PARAMETERS- //) . 

HRITEC6 ,6) 

_JLORKAT„(AX, 


158HCARD 1. -80 COLUMN FREE-FIELD COMMENT CARO 

Z58HCARQ. 2,„~8Q COLUMN FREE-FIELD CCMME NT ..CARO 

358HCARD 3. -NAMELIST/BRGTYPE/ TYF,CAV,PS 


_ii58H 

558H 

J658H ___ 

WRITE (6, 8) 

- FORKAT.dXj.. 

158H 

-25.8JH 


. JJCfL-r.-lLF.CR_ BEARI NG_TY.PE_J) (SEE. ABOVE) 

1 FCR EEARING TYPE i (SEE ABOVE) 
JSLF.CR „B£AR1 NG T_YEE_2_< SEE„ABQVE> 


/IX, 

/IX, 

/IX, 

/IX,. 

/IX, 
I 


CAV - 0 FCR UNC AVITATEO FILM 
1 JLC R__C A VI T AT E 0 JF. IJL JJ 


/IX 

_Z1X 


IF CA V=0 PK=.F. (SEE CARD 7) 
_RS - _ OIL ^SUPPLY. PRESSURE , ..P.SI 


358H 

J» 58H 

75 8HCARD A. -NAMELIST/BE AR ING/ L,R,MU,N 

_85j 8H L r-BEARING _LENGTH,_IN. 

958H R - BEARING RADIUS, IN. 

-JL58H HU -JLUaRICANT VISCJOSIT Y., - MICRO RE.Y NS.. 


/IX 

-/JtX 


/IX 

_Z1X 


/ix 

/AX. 

258H N - RCTCR SPEED (JOURNAL PRECESSION RATE ,RPM) /IX 

„358HCARD-_5«— -NAMELIST/ ECR ATIOA. ES,EF . /AX 

458H ES - INITIAL JOURNAL ECCENTRICITY RATIO, ES>0./1X 

J55&H EF._._.- -FINAL .JOURNAL £CCENTRIC.IT.Y_RAII.O ,__EF<l.iL-/IX 

658HCARD 6. -NAMELIST/CL E ARNC/ C(I),NC /IX 

.75F H Cl JUL -_C LE AR A NO E.. V.A L U E S , _JJu -_0<J < 5 IX X. ; 


NC - NUMBER OF CLEARANCE VALUES 


658H 

WRLTEC6 ,.7Jl 

FORMAT ( 1 > , 

>158HCARDL_7^_-NAHELLS.T/FLGTSEM/_.CS,PC,P.K ,.FP_ 

258H CS - PLOT CONTROL, .T. IF PLOT DESIRED, 

.358 H QJH £R W LSE. ... F._ 


) 


_ZiX,_ 


458H 

_55.8P 

€58H 

_7 5.8 H S.A M PJ. E _Q A T A 


/IX, 
ZIX, 

PC - .T. IF DAMP PLOT CESIREO, OTHERWISE .F. /IX, 

PK - *.T* ...IF._STIFF_. FLOT_ DESIRED, OTHERWISE .F./1X, 

PP - .T. IF PRESS PLOT OESIREO, OTHERWISE .F. /IX, 
/AX $ 


858HCCMMENT CARO i 
_S 5 8.HC CM M E NT. C A R D- JL 


158 H JBRGTYPE T YP = 0 , C A V= 1 , FS=0 . 0 $ 

-25 8H_iBE.AFING. L= Q . 9 Q , R=2 . 55 , KU=Q . 282 , N=i68Q Q$_ 


/IX, 

_/.lX* 


/IX, 

/.IX, 

358H SECRATIO ES-0 • 1 ,EF= 0.9$ /IX, 

A 5 8H_SCL£ ARNC__ C (1 ) = . 003, C (2) = .0Q4,C (3 ) = . 005 , C (4L= . 00 6 , NC = 4$ /IX, 

558H SPLOTSEM CS= . T . ,PC= . T . , PK=. T . , PP= . T .$ ) 

JC 

C WRITE OUT INPUT CATA 

X - - 


WRITE (6, 25) C0MENT1,CCMENT2 

25 FORMAT. ( 1H1,1X, 2 (8A10/1X) ) 

IF(TYP.EQ.O) WRITE (6,27) 

I E_t T Y H . E C . 1 ). . WRITE (6,28 ) 

IF(TYP.EG.2> WRITE (6,29) 

_ - - IF (CAV. EC. 0) WRITE (6,31) 


.27. 

_2_8„ 


IF (CAV. ECvl ) WRITE (6,32) 

FORMAT (/IX, 

158HPLAIN BEARING, NC END SEALS OR OIL SUPPLY GROOVE 
FORMAT LZ1 X , , 


/) 


158HBEARING WITH NO END SEALS BUT WITH OIL SUPPLY GROOVE 
25 FJDRMATX/1X, s 


/) 


SQF00670 
SQFD0680 
SQFD0690 
.. SQFO 07 00 __ 
SQFO 07 It) 

. SQFDG720 _ 
SCFD0730 
.SGFDQ74Q. _ 
SQF00750 
.SQFDC760 ... 
SGFD077 0 
_SQF00780 „ 

SQFD0790 
.SQFD0800 .... 
SQFD0810 

_SQFOG820 

SQFD0830 
. SGFDQ84Q._. 
SQFD0850 

. SQFP0 86Q. 

SQFO 0870 

.5GFDQ880 

SCFD0890 
. SQFP090.0 __ 
SGFD0910 

. S.GFQG9211 

SGF00930 

_SQED09_40_ 

SGFD0950 

SGFOC960. 

SQFD097 0 

. SQFD098 0 

SGFDG990 
SGFD1Q00_ .. 
SQFDiOiO 
. SQFD10 20 __ 
SGFD1030 

_$GF01Q 4Q . 

SQFD1Q50 
SQFD10G0 
SGFD107 0 
SGFD108G __ 
SGF01090 

-SGF.0110 Q 

SCFDlilO 
_SQFD1120 ... 
SGF01130 
_SQFD114G 
SGF01150 

__SQF01i60 

SGF01170 
..SQF0I18Q 
SGFD1190 
SQF01200 __ 
SQF01210 

_SQFP.122.Q 

SQF01230 
_SGFQ124G 
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158H0EARING WITH ENC SEALS AND OIL SUPPLY GROOVE 

_F_Q.fi r.ALUXj UHUNCAVITA7EC FIL M/) 

FORKAT (1X,14HCAVITATE0 FILM/) 

WRITE(6,10) L,R,N,MU,PS,ES,EF,NC . 

F0RMAT(////1X,15HEEARING LENGTH=,F8. 2 , 

128H INCHES BEARING RAOIUS= .F8.2.12H INCHES 

21 X | 15HN= ,F6.1, 

328H RPH MU . F6.3.12H HICROREYNS/ 

41X,15HPS= ,F8.1, 

528H_ _PSJ.. ES= ,F_8.3/_ 

61X»15HEF- ,F8.3, 

_728H NC= j 18// L. 

HRITE(6 ,11) (C(J) ,J=i,NC). 

_L0.R F1U.5.6 Xj. 1 9 HC.LE.FR-A N CE. S - I N CHU// ,5 J SJlX J.F7 ^JA.U/11 


\ 18// L 


OiC«JEl/Z5J BJlXjiEJ j,§./L// /j'.L. 


LE/IN 


IN IIIA.LI ZE_.C0NSI.AM5_ 

_NE.Pr.25. __ 

DELE=(EF-ES)/(NE0-1> 

EI=3 *55/1*13 i 

W=N*0. 10472 



O.U _TE_P_LC CP_F.QR_C_LE.i R ANC E —VALUES. 

ao_jLM_JI=l 4 iiC - 

WRI TE ( 6,12) C(II) 

EO R_MA.IJ.1H i,_l 3 X., AH_£=_ r.FJ. 0 ,_4_,ML IN.* tj_ 

11X,94H EO CO 

__2. PMA X .THETA/ 

31X,97H <DIM) LB-SEC/IN 

_k LB/I N* *2 _ PEG RE ESI 

INNER LOOP FOR ECCENTRICITY RANGE 


00 20 1=1, NEO 

_.KK=mi-l.L* .(.NCr. USJ 1=1. 

EO (KK) = (I-i) ’DELE + ES 

IF (.C.Ay , E.C* 0 >. GOTO 26 

CO ( KK) = (MU*(L**3)*R/(2.*(C(II) **3>)> 

__C0(KK)=CC<KK)*(PI/((1.-EC(KK>**2) )**(3./2.) > 

KO ( KK) = ( 2. * MU * W *R* ( (L/C ( 1 1 ) ) **3 > *E0 ( KK ) ) 

_ K.OJKK) r.KO.LKK ) /.(_( 1 . rrAEQJt K KIA.*? 21* * C 2 ) ). 1 

THETAM(KK>-270.443 - 19 1 . 831*E0 ( KK ) + 218. 223* (EO (KK) **2) 

T HETAM (,KK) = THET AM (KK) .“ 1 14. 803 * (£C.( KK1*_*.3> 

THETA=THETAM(KK)/57. 29577751 

PMAX (KKl=“l . 5* ( (L/C (II) ) **2) *MU*N*EO (KK).*SIN (THETA). 

PMAX (KK) =PM AX (KK) / ( (1 .+ EO (KK) *CQS (THETA) ) **3) 

IFJCAV .EC. 1>_.GOJ.O_AD- 

KO (KK) = 0 . 0 

CO (KK) =KL* ( (L/CJII) ) **3 ) *R?PI 

C0(KK)=CC(KK)/((l.-E0(KK)**2)**(3./2.) ) 

IF (T.YP . NE . 1) . GOTO 2 0 ... 

PMAX (KK) =PMAX( KK) /2« 

CO ( KK) = CO (KK.) /.4 . 

KO (KK) =K0 ( KK) /A • 

CONTINUE. 


SQF01250 

SQFD1260._ 

SQF01270 

_SQF01280 

SQF01290 

SQFD13Q0 

SQFD1310 

SQF0132CL 

SQF01330 

SQF01340 

SQFD1350 
_ SQF01360 
SQF01370 

SQF01380 

SQF01390 

SQFD1400 

SQF 01410 

SQF0142Q 

SGF0143Q 

SQf 01440_ 

SQF01450 

_SQFO1460 

SQF0147Q 

SQF01480 

SQF0149Q 

SQF0150 Q_ 

SQF01510 

SQF01520 

SQFD153Q 

SQFD1540_ 

SQF0155Q 

SQF0156J_ 

SQFD1570 

SQFD1580 

SQFD1590 

SQF01600 

SQF01610 

SQF0162Q 

SCF01630 

SQF01640 

SCFD1650 

SQF01660 

SQF0167 0 

SQF01680. 

SQF0169Q 

.SQF 017 0 0 

SQFO1710 

SQF01720 

SQFO1730 

SQF01740 

SQF017 50 

SQFD1760 

SQF01770 
_ SQF01780 
SOFD1790 

SGFD18U0 

•SGFD1610 
SQF01820 
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300 NN=NEO* NC 

DQ-.3 0_1 11=1 1 , NN , NC 

I=III 

— ™ H R I T E { 6 , 1 3) _ E 0 < I ) , C C ( I) , K 0 ( II, P HA X.il I , T.HELAH ( I >_ 

13 FORMAT (5(5X,F12.3,5X)> 

3 Q CONTINUE _ 

100 CONTINUE 

I El . NQT.a_.CS )_G OJLQ 10 QQ _ . 

NPl=hC*NEO+l 

J4P2 = NC.?_KEG+ NC 

DO AD.O I=NPi,NP2 

... 4 00 EG 111=1. a. . _ . - 

NPi = NPi + NC 

. NP 2= KP. 2 * NC .... 

DO 500 I-NPi,NP2 

— 5.0.fl_EO 1 1 ) = Q-*i66_6Z__ . . _ 


SQF0163Q 

_ SQFD1840 

SQFO1850 
_ _SGFD1860 
SQFO1670 

SQF 01880 

SQF01890 

SaFQ19Q0 

SQF01910 

J5QFD1920 . 

SGFQ1930 

{ SQF01.940 . 

SQFD1950 
. JS.QFD19.63. 


NP=RC*NEO 

—£ ALL-PL CJ.T.ERCNP ,L ,R*MU, N ,PS,QRING., GGMENII,lCQHENI2 ,NC,C_,: 
1E0,C0,K0,PMAX,PC,PK,PP,TYP,CAV) 

-JK -=E 


1000 GOTO 

„393S SIDE. 

END 


90 0 


SQFO1970 

SGF0I98 G_ 

SQF01990 

SQFD20a0- 

SQF02Q10 

S.QFD202.Q... 

SQF02030 

SQF.02fl.4Q.. 

SQFD2050 



1EO,CO,KO,PMAX,PC,PK,PP, TYP,CAV) 


C 

__c 

c 


SUBROUTINE PLOT TER (NP,L ,R , KU , N , PS ORIS G , COHENTi ,CCHEN.T2 , NC ,C, SGFD20 60 

SQFD20 7 0 

SQFD208 0. .... 

SCF02090 

..SGFD21G0 .... 

SQFD211J 

— _SGF_02120„ 

SGF02130 

SQF02i<tQ .... 

SQFD2153 

S.QFO2160 

SQFD2170 

S.GE0218Q .... 

SGFD2190 

SGFD220 0 

SQFD2210 


THIS SUBROUTINE FLCTS CO,KC,ANO PMAX AS FUNCTIONS OF EO 
. E ACH_ S£ T..-OF . CLEARANCE. VALUES . RE GU IRES- CNE BLOCK . 

^DIUENSI C>. C (5) ,£OMENXi( 8) ,C0MENT2 (8X, CLABEL <2)_.___ 

INTEGER COMENT1 ,CONENT2 , CLABEL, TYP 

INTEGER CAV . __ 

REAL L,Ml,N,KO,MUl 
-LO G I C AL ... C RI N G , P C , P K , P R. .... , 


DIME NS I CN £0(135) , CO ( 135) , KO ( 13 5 ) , PMAX (135 3 

JDQ34MON.. JNEQ . _ 

DATA CL ABEL M) , CLAEEL (2 ) / 

Jti 0 HC = K# i D HI LS _ 


_NNF_=NP/.NG 

DELE= ( EC ( NP ) -EO ( 1 ) > / NNP 
_NEH=NNP. 


MUl = MU*t.0E6 
JCflL^QX , -PCJ^GO TJD-.UB1- 


__ploj_e;o._ 


JAXXS=1. 


CALL GRID<L,R»N,HUl,PS,CCMENTi,CCMENT2,JAXIS,EO,CC,KO,PKAX, 

-__1NP., X YP » CAV >_ . , ,,, .... . 


NP i=NC*KEO+i 
.NP2 = KC«NEO+NG. 


DO AGO I=NP1,NP2 

kOJD CO m=o.,.c. 

NPi=NPi+NC 

. NP2=NR2.+JvC. 


DO 500 I=NP1,NP2 
500. .. ...C 0 LI )_= 1 . Q_-. 


C 

_C_ 


00 30 1=1, NC. 

..CALL JLI11EjC£J 0XIUJCXL(XJ^N N.FUflGfJO-,. 01. 


SQF02220 .... 

SQF02230 

SQF022A0 

SQFO/250 

5QFD2Z60 

SGFD2270 

3QF.022 8 0 . _ 

SQFO2290 

-S.GED230.Q 

SQFD2310 

SGF.Q2320 

SQF 02330 

S.QF_023AQ 

SGFD2350 

Jj.G£023.6G _ 

SQFD2370 

SCF0238 0 

SGF02390 

JSQF.O2A0 0. _ 

SGF02410 
SJJE02A20 .... 


C 

X — 


. D£ C I D E_ H HE T.H£R_L A B EL I NG SHOULD BE OVER CUR V£_.0R_ 
AT END OF CURVE. IF MAXIMUM ECCENTRICITY £ 0.7 
_RL ACE— LABEL HG .OVER .CURVE, -QXHERKISE.J LEND 


P EJ.EJCJLNE1 . *2. > G CT G_6 Q_ 


DETERMINE _XH£ X. ANC_Y C CCRCINATES . 0F_IHE. LOWER ... 


C 

_C~ 


LEFT HAND CORNER CF THE LABELING AND THE ANGLE 
I X .MAKE Su. HI I H_XH£_> - A XI S __ 


.XX=CEO ( N£)Y.Q . 16667-)_*-JUjL 
III=NP-NC+I 
..YY= C OCX 1 1) . + fl.i. 


SGFD2430 
._SGE.024AO ._ _ 
SQFDZkSQ 
_SGF02A6Q _. 

SQF02A70 
».SGFD2A.80„ 
SGFO2A90 
_SGFD2500_._ 
SGF02510 
_ SQF02523 .... 
SOFD2530 

„SQFQZ5*iQ 

SQF02550 
SGFD2563 . 


C 

JC_ 


THETA= A T AN( (CO < III 3-CCC III-1G*NC> ) /< (EO (III)-EO <III-i0*NC>> *6. 0) ) SGF02570 

_.THETAi = TH£TA«57.2958 .... .... _ _.$QF02580 

CALL SYMBOL (XX, YY, 0*1 ,CLA BEL , THET Ai , 13 > SGF02590 

SQED260.GL 

DETERMINE THE X ANC Y COORDINATES OF THE LOWER SQFD2610 

-LEEX_ HAND— G QRNER^XF„T.faEL.-CLEARANCE-..VALU£_JX)_BE SQFD2620 
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PLOTTED 

XX=XX*0. 2 6»"C OsTf H E TA ) 

jCY? Y Y+ 0 • 26* SI NJLTJS'E T A ) 

CV=1C00.0*CU) 

..C.ALL__NUKEER_< XX, Y Y ,.0 . 1 ,C Vj THET Al_, 4HF5_. 2 )_ 
GOTO 30 


OETERMINE THE X ANC Y CCCRCINATES OF THE LOWER 

LEFT-HAND. CORNER. OF._THE L.AEELI.NG.ANO _IHE_ ANGLE 

IT HAKES WITH THE X-AXIS 

7x=7 c"( Th^mTneo- iY/ lT+i»' *e . io 

JtHrJ (JNLECrJU LZ\ !.M C + 1± 

YY=C0 (III) + 0.1 

.IP=(IFIX(0, 13129/DELE)) *NC + III . 

THET A=AT AN( (CO (IP) -CO (III ) >/ ( (EO (IP)-EO (III) ) *6 .0 > > 

JLH EXA1 =IH ET A * 5 7, 2 9 5.8 

CALL SYMBOL (XX, YY, 0.1,CLABEL,THETA1,13) 

DETERMINE THE~X A N FT COO R C I N A T ES Tf'Th E LOWER 

..LE£X_HAJs(LiLOJ?.NE^.QF__T_HE__CLE.AR.ANC.E_.VALyE_TO_.SE_ 

PLOOTED 


XX=XX+ G »26*COS (THETA) 

X.Y=JLY._t_JL».2^*.51NiJ li£.T A) ___ 

CV=iOOQ.O*C(I> 

CA LL. Nil H..EER (XX, YY , 0 .1 .CV .THETA 1.AHF5 *_2J 

30 CONTINUE 

CALL__P.L. C I .( 1 5.0 , J ,jQ.,-6 ). 

101 IF ( , NOT • PK) GOTO 201 

J.A XXS= 2. 

CALL GRIC(L,R , N,HUi,PS, CCMENT1, CCMENT2 , JAXIS,EO,CO, KO ,PKAX, 

A N.P. ».I Y.P j .C.A Vi - ' 

NP1=NC* NEO+ 1 

NP.2.=JNC_*HEQ.+ NG ; 

DO 401 I=NP1,NP2 

JtDJ XO.Ulr OjuIT : __ 

NP1=NP1+NC 

.NPj2=KP2+JiC - 

00 501 I=NP1,NP2 

.5.0.1 KO.( I ).= 1 . 0 

DO 40 1=1, NC 

CALI L I LEIE OJD j KO.(I.).,NN P , NC. 0. 0) : 

IF (EO(NP) • GE. 0.7) GOTO 70 

_XX= E C ( NP) *6 . Q._tD. . 1. 

III=NP-NC+I 

YY=KO(.III>..+ 0.1 

THETA»ATAN((KO( III)-KO< III-10*NC) > /( (EO (III ) -EO (III-10*NC) ) *6.0) ) 

T.HE.T.A1=THEI A*57.. 2.950 

CALL SYMBOL ( XX, YY, C . 1 ,CL ABEL, THETAi,13) 

XX=XX+ 0.26*COS (THETA) . 

YY= Y Y + 0.26*SIN(THETA) 

_CV=1000.0*C(I) ... . _._ 

CALL NUHEER (XX, YY, 0.1 ,C V , THET A1 , 4HF5.2) 

.GOTO. 4Q 

70 XX=EC(NC*((NEO-i)/2)+l) *6.0 

III? ((NEC-.i) /2)*NC*I 
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IP=(IFIX(0.i3129/OELE))+tvC + III 

jrx=KOiiii)*o.i , 

THET A= ATAN < CKO ( IP) -KO (I II) ) / ( (EO (IP) -EC CII I > ) *6 . 0) ) ■ 

THETAlr THET A*57. 2958 , 

CALL SYMBOL (XX, YY, 0.1, CL ABEL, THET Ai, 13) 

XX=XX_+J3 . 26* COS < THETA)... 

YY=YY+0.26*SIN( THETA) 

jyr=iooja*jQ*£i.i) 

CALL NUMBER (XX, YY, 0. 1 ,C V » THETAi ,4HF5.2> 

A 0 C ON T I N U £ . . 

CALL PLCTC15. 0,0. 0,-6) 

20.1 IF ( .kO.L*_PPJ -RETURN. 

JAXIS= 3 

CALL .GRin(L,R,N,MUl,P.S»JC.C>!£NIi,CCMEN12.,.JAXIS.,E0,CQ,K0.^KAX., 

1NP , TYP , C AV) 

NP1 = NC + N EO+JL __ * - 

NP2=NC*NE0+ NC 

DO- _4 C 2 _ J « NP 1 NP2 

402 PMAX (I ) =0 *Q 

JNPl^ NP±+Js C ; 

NP2=NP2+NC 

D XL.5 0 2-J = NP.l., NP_2 

502 PMAX (I) =1.0 

D.O„5.0_I.=1, NC : 

CALL LI NE (E C ( I) , PMAX ( I) ,NNF,NC,Q,Q) 

IEILGJLN£T._. GEfc. . (U7J_JG.QT O _8JL ; 

XX=EC (NP ) *6 * 0 +0.1 

* IXI=NP-NC_.+ JE. 

YY=PMAX ( III ) + 0.1 

.THET A=ATARL( FHAX (III) -P MAX NC.) )/ 

1 ( (EO (III )-EO (II I-1C*NC) ) *6.0) > 

XHE.IA1 = JJtiEI A *57. 2.9 5 0 

CALL SYMBOL (XX, YY, 0.1, CL ABEL, THETA1, 13) 

XX=XX+0...26*.COS (THETA) *. 

YY=YY+0® 26* SIN (THETA) 

cv~iooo..o*c ci> 

CALL NUMBER ( XX, YY, G.1,CV, THET A1,4HF5. 2) 

GOXCL5 0 * 

80 XX=EO(NC*((NEO-i)/2)+l) *6.0 

II IhjLC N £ C=1 YJ 2) * N C + 1 

Y Y=PMAX (III ) + 0 * 1 

IP=_(IFXX(Q. 13129/DELE)) *NC +III - 

THETA= ATANC < PMAX < IF > -FM A X (III))/( (EO (IF) -EO (III ) ) +6.0 ) ) 

tHE!A.l= IHEIA* 57.. 29 5 8 — 

CALL SYMBOL (XX, YY, 0.1, CL ABEL, THETAI, 13) 

XX=XX+0.26*COS(THETA) 

YY=YY+G.26*$IN( THETA) 

„ CV=1G0 0 • 0+C (I) . 

CALL NUMEER (XX, YY,0.1,CY,THETA1 ,^HF5.2) 

50 CONTINUE- . ... — 

CALL PL CT( 15.0, 0*0, -6) 

RETURN - . 

END 


SQFD3210 

SQFO3220 _ 

SQFD323 0 

SQFO 3240 .. 

SGFD325Q 

SQF0326Q 

SGFt/3270 

SQFJD3ZSQ 

SGFD329Q 

SQFD330 0 __ 

SQF03310 

SQFD332Q 

SQFO3330 

SQFD334Q. _ 

SQFD3350 

SQED.336D 

SQF03370 

SQF.O3380 

SQF03390 

SQF.0340Q 

SQFD341G 

SQFD342Q. . 

SQFD3430 

_SQF.D344Q 

SQFD345Q 

SQF.Q346Q 

SGF0347G 

S QF 0.348 

SQF0349Q 

SGF035G0 

SQFD3510 

S0.FD352Q 

SQF03530 

SQFD3540 

SQFD3550 

SGFD3560 

SQF03570 

SQF:D3580„__ 

SQF0359 G 

SGFD360Q „ 

SQFO3610 

„ SGFD362Q ... 

SQFD3630 

SGFD364Q 

SGF03650 

SQF0366Q „ 

SGF0367Q 

SQF03680 

SOFD3690 

SQFD37QQ 

SQF0371Q 

S0F0372Q 

SQFD3730 
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SUB.RCUTIN£_GRID(L,R,N,HU1-,FS,COMENT1.,CC('ENT2,JAXIS,EO,CO,KO,PMAX, SQFD374 0 ... 

1NP,TYP,CAV) SQFD3750 

. . .... _ SQF 03760 

THIS SUBROUTINE PLOTS THE CRIO ANO LABELS THE PLOTS SOFO3770 

i SQFD37 80 _ 

01 HE LSI CN C0MENTK8) ,COXENT2{8> , £0(135) ,00(135) ,KC(135) , SCF03790 

1PMAX (135) ,LEGEND1 ( 2 ) , LE GEND2 ( 2) ,LEG£N03 (2) ,LEGEND4(2) ,L£GEN05(3) SQFD360Q_. 

ZHEA0ER1 ( 2) , HEADER2 (3) , XL ABEL (2) , YLABELi (3) , YLABEL2 ( 3) , SQF03810 

3 YL A BEL 3 (3)_. _ ..... SQF 0382 0 ._ 


REAL L , I* , MUi , KO 

INTEGER.. C AV _ , 

INTEGER TYP, HEADER 3(3), HE ADER4 ( 3 ) ,HEA0ER5(2) ,HEA0ER6(2) 

INTEGER-.C0MENT1 »CCI'ENT2 , LEGEND! ,LEGEN02 , LEGENDS, LEGEN04,.. 

1LEGEND5 , HE A DERI ,HE A0ER2 , XL ABEL , YL A6EL1 , YLA8EL2, YLA8EL3 

COMMON . NEO 

DATA LEGENDl(l) ,LEGEND1 (2) ,LEGENC2(1) ,LEGEN02(2) , 

1LEGEND3 (i) ,LEGEND3.(2) ,LEGEN04(1) ,LEGEND4 (2)., LEGENDS (1).,._ 

2LEGEN05 (2) ,L£GEND5 (3) »HEADER1(1) ,HEA0ER1(2) , 

3H£A.DER2.(.l) ,.H£AD£R2. (2.) ,HE ADER2 (3 ) , XLABE.L.(1) , XLABEL (2)., 

4 YLABEL1 ( 1) , YL ABEL1 (2) , YL ABEL 1 (3) , YLABEL2(1) , 

5YLAEEL2.C2L,.YiABEL2i3) ,YLABEL3 (1) , YLABEL3 121 , YLABEL3 (3)/.. . 

610HN= , 10H RPH. ,10HR= I,1QHN. , 

7J..0»L= 1,1 OHN.. ,10HPSS ,10H_ ,PSI. 

810HMU= , 1 OH MCROREYN, 10HS ,10HSQUEEZE FI, 

910 HLJl...DLAi!PER ,10 HNC. .SEALS. ..C,10HR_OIL..SUFP,10HLY__ GROOVE,. 

110HECCEKTRICI , 10HTY - (DIM) ,10H CAMPING, ,10HCO- (L8-S, 

210BEC/IN) ,1 OH STIFFNESS , ,10H . KO ..r._lLE/, l.OHIN) ,. 

310HHAXIHJH PR,10HESSURE - (,10HPSI) / 

U.AT.A...HEACER3(1) .HEACER3 (2) , HEADERS (3) ,HEA0ER4(1) , 

1HEADER4 ( 2 ) , HEA0ER4 (3) , HE A0ER5 (1 ), HEADERS (2) , 

2HEADER6.(1) , HEAD.ER6.(2)./ 

310HCIL SUPPLY, 1 OH GROOVE -,10H NO SEALS, 10HGIL - SLIP PL , 

410KY— GRCOV E_ ,1. OH AND • SEALS ,10H. RO CAV ,1 0 HI T A.TION ,.. 

510H CAVIT ,10HATION / 

IC Y OLE S = fl_ " 

LOEXP=0 

G Q_ID_U..O 0 , 2 0 0 , 3 0 0 ). , J AX IS.. 

100 CALL L0GSCAL(C0,NP,1,8. 0 ,LCEXP, I CYCLES , NP*2) 

u_G.0TD.JL5O 

200 CALL LOGSCAL (KO ,NP,1 , 8» 0 , LCEXF, IC YCLES , NP+2) 

.._.G0I0_15 0 L_. 

300 CALL LOGSCAL (PHAX , NF ,1, 8 . 0 ,LO£XP , ICYCLES,NP+2) 

IRQ CQN.T INUE I 

NNC=(NP/NE0)+1 


SQF03830 

SQF.03840 _ 

SQF03850 

S.QED.3 8 6 0_ 

SQF03870 

SQFD3880_ 

SQF03890 

SQED390Q. 

SQFD3910 

....SQF.D3920L 

SQF03930 

S QF.0 3.9.4 0 . 

SQFD3950 

SQF 03960.. 

SQF03970 

SQFD3980.. 

SQF03990 

SQF Q4 000.. 

SQFD4010 

SQF_D402(L. 

SQFD4030 

SQFD.4.04.0.. 

SQFD4050 

S.QFD4060.. 

SQFD4070 

SQFD408Q 

SQF04090 

SQFD4100 

SQF04110 

SQF04120.. 

SQFD4130 

SQFD4140 

SQF04150 

... SQF0416Q 

SQF04170 


CALL. AXISiiO. 0,0.0 , XLABEL, -20, 6. 0,0.0 ,£D(NR+1) »E0 (NP+NNC) ,16.667) SQFD418Q ... 

GO TO (1G1 ,2 01, 301) , JAXIS SQFO4190 

.101 CALL LOG AXIS (0. 0,0.0, YLAEEL1, 26, + 8. 0,9 0,0, L.OEXP.ICYCLES) ..SQF0420 0 ... 

GOTO 160 SQFD4210 

.201 CALI L 0 G A X I S (.0 • fl ,0-.J ^YLA EEL2 , 23, + 8 . Q. , 9 Q . 0 , L 0 EXP ,1 CYCLES). ..SQFD4220 

GOTO 160 SQF0423Q 

-.301—- CALL LOG AXIS (0. 0,0.0, YLAEEL3, 24, + 8. 0,9 0.0, LOEXP,ICYCLES) SQFD4240 ... 

160 CALL PLCT (0.0, 8. 0,3) SQF04250 

CALL i>LCT(6. 0,8. 0,2) SQFD4260... 

CALL LOGAXIS(6.0,0.0,2H , Q»+8. 0 , -90. 0 jLOEXP , ICYCLES) SQF04270 


PLOT LEGENDS AND HEADINGS 
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CALL SYMBOL (1.29*8.96,0 .2 1 » HEADER1 , 0. 0,19) 

IF(TYP.EG.O) CALL S YMBOL ( 0 . 30 , 8 . 67 . 0 . 2 1 .HEA0ER2 , 0 . 0 , 30) 

SGF0431Q 

SQFD4320 


IF (TYP.tG.l) CALL SYMBOL ( 0 .30 , 8 .57 , 0 . 21 ,H£ADER3 , 0 . 0 ,30 ) 
IF (TYP.EC.2). CALL. SYMBOL <0 . 30 , 8 .57, 0 . 21 ,HE A0ER4, 0 . 0 ,3 0) _ 

SQF0433Q 

SGFD4340 


IF(CAV.EC.O) CALL SYMBOL ( 1 . 20 ,8 . 18 , 0. 2 1 ,HEA0ER5 , 0. 0 , 20) 
IF(CAV.EG.l) CALL S YMBOL ( 1 . 20 , 8 . 18 , 0 . 21 ,HEAOER6 , 0 . 0 ,20) . . 

SQFD4350 

SQFO4360 


CALL SYMBOL (4. 0,1. 0,0. 10 , LEGEND 1 ,0.0,16) 

SQFO4370 


CALL SYMBOL (4.0 .0.8. 0.10, LEGEN02, 0.0, 12) . ....... 

SQFD4380 


CALL SYMBOL (4 . 0 , 0 . 6 , 0 .1 0 , LEGEND 3 , 0. 0 , 12 ) 
CALL SYMBOL <4.0 ,0.4.0.10 , LEGE NO 4 ,0.0,16) 

SQFO4390 

SCF04400 


CALL SYMBOL (4.0 , 0 .2 , 0.10 , LEGENDS , 0 .0 , 21) 
CALL NUMEER <4 . 26 . 1 . 0 . 0. 1 0 . N . 0 . 0 ,4HF8. 1 > 

SQFD4410 

SQFO4420 


CALL NUMBER (4.26,0.6,0. 10, R, 0.0, 4HF5. 2) 

CALL NUMBER <4. 26. 0.6.0. 10.L.0.O .4HF5.2) .... . 

SQFD4430 

SQFD4443 


CALL NUMBER(4. 26, 0.4,0. 10 , PS, 0 . 0 , 4HF7 ,4) 
CALL NUMBER ( 4 . 26.0.2 .0. 10.MU1.0.0.4HF6.3) 

SQFD4450 

SQF04460 


CALL SYMBOL (-0.43, -0.75, 0.10, COMENtl, 0.6, 80) 
CALL SYMBOL (-0. 43 0. 90 . 0 . 10 , C0KENT2 ,0 . 0 .8 0 ) 

SGFD447 0 
SGF0448Q 


RETURN 
END 

SQF0449D 
SCFD450 0 



S UB R.CUX I JL 0 use AJL .1 Jf , K»_K , S., M3 * QJf SLL E S j. K N 1 SQFD451Q 

C ‘ SQFD4520 

£ THIS. SUBROUTINE TAKES .AN. ARRAY_.OF . DATA *T*_AND_SCALE.S_ JIS..VALUES F0SCFO453Q 

C A LOGARITHMIC PLOT. THERE ARE TWO ENTRIS SQF0454Q 

& L0-6.SC.AL_I S._U.SEC._If_ J_HE_.RANG E _OF OATA_ VALUES, ARE UNKNOWN _ _SQFD4550 _ 

C SQFD456Q 

C. LOj».C^JL_I^_y.S£0-J.F_IJE_RA.N6E--0.F_-0.AlA_..A8£.j<NO WN_A NO_C.AN_.BE_.PRO V ISQFD4570. _ 

C SQFD45BQ 

C .INPUT -PARAMETERS FOR 6CTH RCUTIN SQFD4590 

C Y= THE ARRAY IN WHICH THE DATA TO BE SCALED IS STOREO ~ ' ”SQFO460 «) 

C K=.TJHE_ NUMBER . OF. JHE POINTS IN THE ARRAY SQFD4610 

C S= THE LENGTH OF THE AXIS IN INCHES SQFD4620 

C K=A REPETITION FAQTCR__FOR_ACCESSING DATA IN _Y_ SQFD4630 _ 

C INPUT PARAMETERS FOR LCGCYC BUT OUTPUT PARAMETERS OF LOGSCAL ' s’cFD4640 

C_ CYCLE S_= T FE. NUMBER .CP. CYCLES. OF REPETITIONS OF THE GRAPH SCFD465Q 

C LO=THE LOWEST EXPONENT OF OF THE CATA SQFD4660 


C IE THE DIFFERENCE BETWEEN HIGHEST AND LOWEST EXPONENTS + ONE SQFD468Q 

__C I E. I F Y (I ) IS GREA TER THAN 1.0*10 »»4 THEN LO =4 SQF04690_ 

C INTERNAL VARIABLES SQFD4700 

C EI.-S.IHE_ M AXI HUN _VAL LE_CP_T HE.. DATA SQF 047 1 0 _ 

C AHI AND ALO ARE REAL VALUES FOR INTEGER HI AND LO SQFD4720 

CL CL._5.A_ SCALE -FACTOR SQFD4730_ 

C ♦♦»**»»******»»»»*»»****4*»»»**»»#*»*»*»»*****»*»*»»»**i»*»*M *****SCF04740 

INIEGER HI 4 ...CYCLES __SQFD4750_ 

DIMENSION Y <NN) SCFD4760 

AH I =-J. ,Q.._L*3.21 SQF 04770 

ALO=10.**32i SQFD4780 

C SQFD4790 

C FIND THE MAXIMUM OF THE DATA SQFD480Q 

C_, SQFD4810 

5 DO 10 1=1, N,K SQFD4820 

AHI = AHA X i.< AHI Y ( I.)_.)_ SQF 0 483 0 _ 

ALO= AMI N1 (ALC , Y ( I) ) SQFD48‘40 

10 CONTINUE . SQFD4850 

C SQF 04 860 

£L .C.OHEULLS—L C HI . AND C YQJLES ; SQFD4870 

C SQFD4880 . 

.A.LO=AL.PG10 < ALO) SQFD489Q 

AHI=AL OGiO ( AHI) SQFD4900 

LOsIF.IX.CALO) S0FD4910 

IFTALO.LT. 0.0) L0=LC-1 SQFD4920 

: HJsIEIXJA HI. Li 1 SQF 0.4 9 30 

CYCLES=HI-L0 SCFD4940 

_JGQ.lJO.JU ; .SQFD495Q 

‘ C SOFD4960 

C IH1S... ENTRY. POINT.. IS. USED. .WHEN. LO... AND. C.Y_CL.ESl.ARE„.KNOWN SCFD4970 

C SQF04980 

ENTRY. LCGCYC. SQFD4990 

11 CL=S/CYCLES SQFD5000 

00.-2A 1=1, . N,..K SQFD5010 

C SGFD502Q 

C _GET_ LOG 10_OF THE DATA AND .AJ3JUST TO SCALE _ - _ SQF05030 

C SQF05040 

IF_C ..YU) .LE.O.) PRINT. 30 SOFO505Q_ 

30 FORMAT ( IX, ’NEGATIVE OR 0 VALUES HAVE BEEN ACCESSED BY THE SCALING SCF05060 
_iSUBRQUTINE**.*f***. ,, *. 1 L*.-THE_ ABSOLUTE VALUE.OF..THE.VALUE .HAS BEEN USSCFD5070 
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iED FOR PLOT TING PURFOSES) 

YXl ) fAL OG i 0 < AES CY-II.) ) )— 

20 Y ( I ) = <Y (I ) *L C) 4 CL 

- XIN+1) = C,0 

Y (N + 2) = i • 0 

RET URN ; 

END 


SQF05080 

SCFD5Q90 , 

SQFD5100 

_SGF 05110 

SGFO5120 

SCFO5130 

SGF0514Q 



— SU BROU TINE. -LOG AXIS (Xp_Y..» IBCD,.. NCH ,. SZ , -THT A LOEXP, JCYCLES) SQFD515Q. 

THIS ROUTINE CREATES A L CG ARITHNIC AXIS WHICH MAY EE HORIZCNTAL SQF05160 

(ALCNG THE LENGHTH CF . THE .P AFER OR THTA = 0.0) OR . VERTICAL SQF05170 

{ACROSS THE WIDTH OF THE PAPER OR THTA =90.0) , OR AT ANY SQFD5160 

OTHEk ANGLE*- ELR.IHER_.IFE. AXIS_CAN_ BE. DRAHN. _HITH_J3ESIRED LABELISQFO5190 

PRINTED SQF05200 

.VARIABLES, f* 1 ^** .***•* f ****** SQF.D521Q 

INPUT PARAMETERS SQFD5220 

_X, Y-.THE. COORDINATES OF . THE BEGINNING POINT.. OF. T.HE. AXIS_ SQFD5230 

I8CD = AN ARRAY OF HOLLRITH CHARACTERS PRINTED AS LABELING SQF052AQ 

..NCHARrNUKBER OF ..CHARACTERS IN IBCD. T.O..FLOT SQFD5250. 

SIZE = THE LENGTH OF THE AXIS IN INCHES " SQFD5260 

— THE I.A_- . _1HE . . A NG.L E- .A T_ .WH I C H...THE...AX IS ..I S_DR A W N SQFD527 0. 

LOEXP-MI NIMUM EXPONENT{ POWER OF TEN) CF ANY VALUE TO BE PLOTTED SCF05280 

— ICYCLES.- NUMBER. OF . .REP ITJJLICNS OF. -THE AXIS- J.Q_B£._PLQTT.EQ SQFD5290. 

INTERNAL VARIABLES SQFD530Q 

—F.INE. DETERMINES FINE .OR COARSE -TICK. MARKS SQFD5310.. 

HGT-CONTAINS THE HEIGHTS CF THE TICK MARKS SQFD5320 

— COQRR_.r.£CNT AINS...LG£S_US£C...JO_MO.VJE__J.HE_JBEN _T.H£_C.08RE.GI_L£NGHI SQFJD533Q- 

WHILE THE GRAFH IS BEING DRAWN SQFD5390 

»»» * ** » ♦»»» »«»» •* * * ***+j?J£+* +**+*:* *** ?„****- SQF D 533 Q_ 


DIKENSI Oh COORD (65) » HGT(65>, IBC0(6) 
..LOGIC 


: KOJE_XE KJEEQ_TH£_AilSOLUTE_VALU£S..JDE__J5LZ£*_TJi£JA|„NCHAR__ 


„SIZE=SZ 

NCHAR=NCH 

Til £ T A* T h JA 

2=COORDU) 

-2 1 s. AL 0 G. C£LQ All 11 1/1 0 « 
IF ( Z «EQ« ZZ) GO TO i 


C THIS PROCEDURE FILLS THE ARRAYS KITH NEEDED INFORMATION 

c ; 

ITEN=6H 1X10 

K*1G ; 

INCR=i 

DO .3 X=1*.6SL __ 

K=K+INCR 

A=RLCAT 00/10. : 

COORC(I)=ALOGiO(A) 

=l*.04 

J=K-(K/5>*5 

IF ( J..EQ • Q ) HGT. II . 06 _ 

«J=K“ (K/ 1 0) * IQ 

IF LJj,EQ. 0) . HGT II)-* 100 

IF < K .EQ * 50 ) I N C R - 2 

5 — CONTINUE : 

C 

C THIS -SECTION SETS- ThE -VA LUES FOR -THE -DIFFERENTJPPIIONS _IN_.TICK . 

C HARKS LABELING ETC* 


DIRECT = 1 • 

IE1NCHAR.LTL..0J DIRECT^- 1.) 

TICK ^DIRECT 

JEIIHEJ.A..LT *U .l _IICK_s{ -T.XCKi__ 
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FINE=. FALSE • SQF05720 

IF_(SIZE,LX,.Q . ) FIM_J=JlECE.*___ S.QF D5730..._ 

THETA=AES(THETA) SQF05740 

NCHAR_ = IAB$ (NCHAR) _ SQFD5750 

SIZE =ABS{SIZE> SQFD5760 

CS T HETA=COS (THETA* .017455 ) SQFD5770 _ 

SNTHET A = SIN (THET A* . 0174 5 5 ) SCF05780 

_0_ . • _ _SQF05790 _ 

C THIS DRAWS THE AXIS AT THE FROPER ANGLE AND LENGTH SQFD5800 

C_ IT USES_ COORD ARRAY FOR GIVING THE RIGHT LENGTH FOR TH AXIS AND TICSQFD5810 _ 


_ C DIRECTION. QF_TH.E LABELING ANDOF_jHENUMB_ER O F T ICK HARK S _SQFD5830_„ 

C SQF05840 

ClC.LES.Zr SIIE/FL OAT.dC YCLES) SQFD5850_ 

V=. 10* (-TICK) *SNTHETA+X SQFO5860 

. H=TICK* . 1 0*CSTH£T A+ V S0FD5870 _ 

CALL PL CT ( V » H, 3) SQF05880 

CALL. PLOT (X, Y,_2) ~ . SGF05893 

K=5 SQF0590Q 

If ..( F INEIKEJL . SQF 0.59 1 0_ 

00 20 1=1 jICYCLES SQF05920 

■l= fl : SQF 0593 3 

CALL WHERE (C, E, IOUHHY) SQFO5940 

ID J= J+K . SQFD5950 

A=COORD ( J) *CST HETA*C YCLESZ+C SQF0596Q 

B= C 0 OR D ( JA1S.NT HE Lfi* CYCLES! + 0 S QF 0 5 97 0._ 

CALL PLOT (A, B, 2) SQF05980 

V=HGT ( J). * (.-TICK) *SNTHETA + A SQFO5990 _ 

H=TICK*HGT(J)*CSTHETA+B SQFD6000 

CALL PLOT (V,W,1) SQFO6013 

CALL PLOT (A, B, i) SOFD6020 

IE ( J •,L_T..»._6 5.)_G.0_..T.C _10 S0FO603.0_. 

20 CONTINUE SQF0604Q 

IF (NCH ARjiECL* 0) _G0_T_C__82J1 _SCF06050 

C SQF0606Q 

C .THIS. ROUTINE PROVIDES LABELING IN THE PROPER DIRECTION, SCFD6070 

C FOR THE IDENT IFICATICN OF THE TICK HARKS SQFD6080 

C : SQF0609Q 

IQ* 5 SQF06100. 

IUL=_40 _ _ SOFD6110 

310 BB= • 07 * TICK +.21*DIRECT -.07 - ( TICK+OIRECT ) *. 03 SQF06120 

OX=CSTHETA* (-.2 8) - E£*SN THETA SCF06130... 

DY=SNTHETA* (-.28) +E8*CS THETA SQF06140 

XA=SIZE*CS.T HEIA + X+EX S0F0615 0 ... 

YA=SIZE *SNT HET A + Y + CY SQFO6160 

CALL SYHBOL ( XA, .YA, .,14., .I.TEN,_THETA,_6J : SQF06170 

11= ICYCLES+LOEXP SQF06180 

XA=CSTHETA*_. 72-SNTHETA* . 05 + XA. CQFO6190 . 

YA=SNTHETA* . 72FCSTHETA* . 05+Y A SQF06200 

P0WER?.2HI.3 : SQF.D6210 .. 

CALL NUMBER (XA, YA, .10, II, THETA, POKER) SQF06220 

Ilr.ICYCLES— SQF 0 £2 30 ... 

C SQFD6240 

c. THIS MOVES ..THE PEN CONN T H E_ A XI S._P UT T I N.G..T.H E. TICK MARK .LABELING .. SQFD6250 

C AT THE CORRECT POSTION SQF06260 

c. - SQ.F06270 

320 K=IO SQF06280 

11= 1 1-1 ..SQF06290 
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INCRr-5 

i=xui 

QA = (CO ORO ( I ) + FLOAT ( 1 1 ) ) *C YCLESZ + .25 

X A = Q A* C S T KET A+ X + 0 X 

YA==QA*SNTHETA+Y + OY 
K = K - i 

^ c i # £ qT^oT *~I N Vr 

I = I.+ 1 N C R 

IF (I.GT.O) GO TC33G 

GA = F L 0 AT (II) ♦CYCLES 2 _ . 

X A = Q A*C STH E T A + X + DX 

JTA _ = C A * S NT H E T A + Y + OY 

CALL SYMBOL (XA, YA, .14, IT EN , THETA, 6) 

Il=il+L OEXP 

XA=C STHET A* • 72-SNTHET A*.Q 5+ XA 

_Y_A=SK I H E T A * . 7 2 + C 5 IH E I A ♦. 05+ Y A 

POWER =2 H 13 

— CALI — NUMBER! X A YA, — *10., IL.,_t_HET A.,_ POWER! 

IFMI.GT.Q) GO TO320 


THIS PRINTS THE DESIRED ALPHA NUMERIC INFCRHATION ABOVE OR BELOW 
_IHE_.AXI.S_AS-.IS .N£EO£D.__TO I CENTIF_Y_ THE_ .AXIS 


__B A=. 5+S I ZE- FLOAT (NCHAR) ... f . 06 

BB -TICK*. 05 +D I RECT * *4 G- • 07 

—X A=BA*ES-T.H E TA-B S* S N I H ET. A + X 

YA =BA*SMHETA + B8*CSTHET A + Y 

_ C A L L . .S Y K£ OL LX A , _ Y A ,_._14. , ...I BCD, THETA, NCHAR). 

RETURN 

_END_ . 


SQFD6300 

SQFD631Q_ 

SQFO6320 

.. SQFD6330 . 

SQF06340 
_ 5QFD635Q_ 
SQFO6360 
___ _SQED637 0 
SQF 06380 

SGFD639G 

SQF 0640 0 

S Q£ D 64 1 0_ . 

SQF06420 

SQFD6430. 

SQF 06440 
__SGFD645a_ 
SGFD6460 
„SQF5647Q_ 
SQF06480 

S.GF_D649.0_ 

SGF06500 

SQFD651Q 

SQFD6520 

SQF06530- 

SQFD6540 

SQF.D655Q„ 

SQF 06560 

SGFD6570- 

SQF0658Q 
SQFD6590_ 


r ~**r»Tc cefzfti lth: iWETT^ 1 ^ 

“ mis prcgram analyzes nrTrrFnyESsy d'atfing and ffetstre 

CHARACTERISTICS OF THE S0DEE2E FILH CAFFES BEARING. THREE 

“Tea r I’NTWK'FTCiirA tt cns - p~ay" _ ee a n aitz ett- 

0 - PLAIN BEARING WITHOUT END LEAKAGE SEALS OR CIRCLH- 

FERE KTT A L~0 ILTUFF CY — TfFO CVE — — — — 

1 - BEARING WITHOUT END LEAKAGE SEALS BUT WITH CIRCLH- 

FERENTIAL OIL SUPPLY GROOVE 

2 - BEARING WITH BOTH END LEAKAGE SEALS ANO S1RCUHFER- 

nrrixirTnirscFPCY groove 

IN ADDITION, THE FILM HAY EE ASSUHEO TC EE EITHER CAVI- 

tAtEO cfi uNTAvrnrrEcr'iF cAvrirsTErTFr tilh is asscket 

TO EXTEND FRCN THETA=PI/2 TC TH£TA*3FI/2, WHERE THETA 

"IS FE'SsCFTC' "FRCF ThE LI NF _ CF"C'E'KTE'R'S IT T RET D IRE C T ITNTF 

JOURNAL PRECESSION, THE EVALUATICN OF THE BEARING CHARAC- 

~ TErTsTITTS - 'A'STiJFET THAT" THE - j OO'RN A C“FREXTSSTS _ S YR C H R CN CTT- 

LY ABOUT THE BEARING CENTER. 

~THE'ToLLWING IS A CESCfilPTICN CF THE TNFll PARAMETERS- 


CARO 1. -80 COLUMN FREE-FIELO CCHHENT CARD 

“TARO" 2T*^80 — CoCuMfr-FR'E^FIELO CCHHENT CARC 

CARD 3. -NAHEL IST/ERGT YFE/ T YP ,CA V ,PS 

TTF - O' FOR E E ATRTN G TYPE' 0 TSEfc ABOVE) — 

1 FOR BEARING T YFE 1 (SEE AeOVE) 

Z FTTR - EETRINTj TYPE 2 TSfct'" ABOVE ) 

CAV - 0 FOR LNCAVITATF.C FILH 

rrCfR" C'AVTT ATED ' FI LH . 

IF CA V=0 PK=.F. (SEE CARD 7) 

PS - TI tT'S'L F PO PRE SSCKE , _ FST — 

CARD 4. -NAHELIST/E EARING/ L,R,HU,N 

C "BEARING LENGTH, IK. 

R - EEARING RADIUS, IN. 

HIT - Tl ERI C ANT" VTSCTCTITYT - FITRCREY NS 

N - RCTOR SPEED (J64SE13 7IEC5S29C5 IlCE.RFM) 

“TAR D — 5 . "-KAHECI ST7ETRATI137~ES7ET — 

ES -.INITIAL JOURNAL ECCENTRICITY RATIC, ES>3. 
ZF~ - FINAL" JOURNAL ellEn (K ill i y kaIIc, tH<i.a 
CARO 6. -NAHEL IST/CLE AFNC/ CTI),NC 

cm- clear atct taoj e st it; — tr<iT5 — — ~ 

NC - NUHBER OF CLEARANCE VALLES 
"CARD 7. -NAkELIST/FLOTSEM/ CS,FC,FK,HH — — 

CS - FLCT CCNTRCL, .T. IF FLCT DESIREC, 

5 T FE ' f T H IS E — TF^ — — — — 

PC - .T. IF DAMP PLOT DESIRED, OTHERWISE .F. 
PK" - TT. IF STIFF.” FLCT CETT REC, OlHfcRWibt ,K. 
PF - .T. IF PRESS FLCT DESIREC, OTHERWISE .F. 

— SATfPLrTATA ~ — — — — — — — - 

COMHENT CARD i 

COHMENT CARC 2 ~ * — — 

SBRGTYFE TYF= 0, CA V=1 , FS-0 . 0$ 

fEfAR I N G L= C-. SO ,’R = 2 . S ? , KU= 0T3 82 , N = l68"TS'01 

SECRATIC ES=0.1 ,EF=0. 9? 

iCLE ARNCTU' ) = . £T0"3 ,C < 2 ) = . 0 04 , C (3 ) s . qO S , C ( 4 ) = . (TC 6 , NU= 4 % 
SPLOTSEH CS- . T . ,PC = .T. , PK~ .T . , FF = «T .S 
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<Oin L6-SEC/XN LE/IK L3/IK**2 DEGREES 
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